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ABSTRACT 


In  the  procurement  of  electronic  systems  by  the  Department  of  Defense, 
electromagnetic  compatibility  has  generally  been  considered  in  the  design 
phase  by  the  application  of  somewhat  rigid  and  rather  general  specifications 
to  the  individual  equipments/subsystems. 

This  report  examines  the  pertinent  stardards  and  reviews  existing  efforts 
to  predict  intrasystem  compatibility  via  computer  aided  techniques. 
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I.  INTRODUCTION 


In  the  procurement  of  electronic  systems  by  the  Department  of  Defense, 
EMC  has  generally  been  considered  in  the  design  phase  by  the  application 
of  somewhat  rigid  and  rather  general  specifications  to  the  individual  equip¬ 
ments/subsystems.  For  example,  MIL-STD-461  places  limits  on  the  allowable 
emissions  from  a  subsystem  (a  "black  box")  and  also  places  limits  On  the 
susceptibility  of  jubsystems  to  emissions  from  other  subsystems.  This 
MIL-STP  basically  attempts  to  insure  a  "quiet  system,"  i.e.,  the  specifi¬ 
cation  controls  undesired  or  noise  type  signals  within  the  system. 

The  problem  here  is  that  the  MIL-STD‘s  (e.g.,  MIL-STD-461,  M1L-STD-826, 
MIL-STD-6181)  are  not  written  individually  for  each  system  under  considera¬ 
tion.  They  are  a  set  of  specifications  which  are  to  be  applied  to  all  sub¬ 
systems/equipments  the  Air  Force  (as  well  as  all  departments  of  DOD  for  the 
current  MIL-STD-461)  purchases.  The  present  subsystem  EMC  specification, 
MIL-STD-461,  seems  to  have  evolved  over  the  past  30  years  and  it  is  dif¬ 
ficult  to  surface  any  rationale  for  the  limits  in  these  specifications.  In 
fact,  it  seems  that  it  is  fundamentally  impossible  to  attach  any  rationale 
to  these  general  specifications  at  all.  To  do  so,  one  would  have  to  assume 
that  one  set  of  specifications  could  insure  compatibility  for  such  a  wide 
t'ass  of  applications. 

To  illustrate  the  deficiencies  in  the  present  approach,  consider  the 
following.  Th **  electromagnet  c  energy  emitted  by  subsystems  will  be  coupled 
to  other  subsystems  through  various  transfer  paths.  For  example,  signals 
on  a  power  line  which  are  not  at  the  power  frequency  may  be  undesirable 
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since  they  may  cause  interference  by  coupling  to  other  subsystems  as  a  result 
of  direct  conduction  through  the  common  power  supply,  radiation  to  nearby 
sensitive  receiving  antennas  or  radiation  to  other  wires  and  equipment 
cases.  Also  the  signal  on  the  line  at  the  power  frequency  will  couple  to 
ocher  subsystems,  possibly  causing  interference.  What  would  be  very  de¬ 
sirable  to  know  are  1)  the  levels  of  the  desired  as  well  as  undesired  emis¬ 
sions,  2)  the  portions  of  these  emissions  which  are  incident  upon  each  sus¬ 
ceptible  device  due  to  transferral  through  each  coupling  path  and  3)  the 
quantitative  susceptibility  of  devices  to  the  incident  energy.  If  one  can 
determine  these  items  then  a  meaningful  determination  of  compatibility  can 
be  made. 

Herein  lies  the  primary  deficiency  in  the  present  approach  to  EMC  of 
applying  the  same  LMC  subsystem  specifications  to  many  different  systems. 
Simply  placing  limits  on  allowable  emissions  and  susceptibilities  of  sub¬ 
systems  will  not  necessarily  insure  compatible  operation  of  these  subsystems. 
In  certain  cases,  the  actual  emissions  from  the  subsystems  (even  though  they 
do  net  exceed  the  limits  set  by  the  MIL-STD's)  when  transferred  through  the 
various  coupling  paths  may  well  exceed  the  prescribed  susceptibility  limit 
(set  by  the  MIL-STD's)  of  a  device.  On  the  other  hand,  the  limits  may  be 
too  stringent  requiring  unnecessary  suppression  Measures.  Therefore,  to 
achieve  compatibility  for  a  system  and  optimize  the  required  interference 
suppression  measures,  one  should  determine  the  energy  coupled  between  the 
subsystems. 

Clearly  a  subsystem  will  be  affected  to  some  degree  by  all  sources  of 
energy  (e.g.,  desired  and  undesired,  wires,  cares  and  antennas,  etc.).  A 
fundamental  requirement  tor  insuring  compatibility  of  a  system  then  is  the 


consideration  of  all  paths  of  electromagnetic  energy  transfer.  Not  only 
must  one  consider  the  totality  of  these  paths  but  the  total  amount  of 
electromagnetic  energy  transferred  must  be  determined;  at  least  the  "sig¬ 
nificant"  energy  must  be  determined.  Obviously,  the  terms  "all  paths" 
and  "total  amount  of  transferred  electromagnetic  energy"  are  only  concep¬ 
tual  topics.  In  a  practical  application,  achieving  these  objectives  would 
He  an  impossible  task.  We  are  only  emphasizing  the  obvious  point  that  all 
sources  of  energy  will  tc>  some  degree  affect  each  subsystem  in  an  actual 
deployment.  This  comment  points  out  that  not  only  should  one  consider 
coupling  of  energy  generated  within  the  system  but  also  energy  generated 
outside  the  system  which  is  incident  upun  the  subsystems  should  not  be  dis¬ 
regarded.  Also  one  must  consider  desired  or  functional  signals  as  well  as 
undesired  or  nonfunction il  signals.  For  example,  one  cannot  neglect  the 
power  supply  fundamental  frequency  (e.g.,  60HH,  400H2)  in  determining  sub¬ 
system  compatibility  since  the  subsystems  will  be  affected  to  some  degree 
in  an  actual  deployment  by  these  signals. 

A  system  can  be  defined  in  many  ways  and  is  subject  to  a  great  deal  of 
variability  depending  upon  the  particular  needs.  One  muy  define  a  resistor  as 
a  "system.''  On  the  other  hand,  one  may  define  the  entire  universe  as  a 
"system."  Between  these  two  extremes  lies  a  compromise  definition  of  a 
system.  For  example  one  may  define  an  entire  precision  approach  radar 
complex  as  a  system.  Within  this  system  may  be  included  the  facilities 
:  the  PAR  site  to  communicate  with  the  outside'  world  such  as  the  VHF 
voice  comaunications. 

Clearly  the  definition  of  a  system  depends  upon  tne  situstion  at  hand. 
Within  a  system  we  'aaj'  isolate  subsystems.  The  definition  or  delineation 
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of  these  subsystems  generally  is  a  compromise  of  two  considerations.  If 
one  breaks  the  system  into  many  subsystems  which  are  physically  small,  then 
the  analysis  may  become  too  unwieldy.  For  example,  one  could  define  the 
subsystems  for  the  PAR  complex  at  the  component  level,  i.e.,  resistors,  in¬ 
ductors,  etc.  However,  if  one  defines  the  subsystems  on  a  more  macroscopic 
"black  box"  level,  e.g.,  transmitters,  power  supplies,  etc.,  then  the  result¬ 
ing  "model"  of  the  system  consisting  of  an  "interconnection"  of  these  black 
boxes  may  not  admit  an  acceptable  degree  of  accuracy  in  the  system  analysis. 
For  example,  one  may  model  a  communications  receiver  as  a  frequency  dependent 
sensitivity  curve.  Yet  this  macroscopic  model  does  not  predict  other  effects 
such  as  intermodulation,  crossmodulation,  etc.  Alternately,  the  receiver 
could  be  modeled  at  the  component  level  to  include  all  these  effects  with 
the  resulting  complex  analysis.  One  could  conceive  of  an  electronic  system 
as  a  set  of  "black  boxes"  which  are  connected  by  wires  (cables)  and  antennas. 
In  other  words,  these  boxes  would  communicate  or  interact  with  each  other 
through  the  transfer  of  electromagnetic  energy  via  coupling  between  wires, 
antennas  and  the  boxes  themselves. 

In  order  to  determine  quantitatively  the  coupling  of  electromagnetic 
energy  between  the  various  subsystems,  a  mathematical  system  model  should 
be  determined.  This  model  would  consist  of  analytical  models  for  genera¬ 
tors  and  receptors  (or  processors)  of  electromagnetic  energy  along  with 
analytical  models  of  the  various  transfer  paths.  The  use  of  the  term  sys¬ 
tem  model  does  not  imply  that  there  will  be  one  model  for  all  systems.  Ob¬ 
viously  this  cannot  be  the  case.  A  system  model  is  merely  a  mathematical 
interconnection  of  the  various  subsystem  mathematical  models  to  represent 
the  actual  system. 
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In  a  realistic  situation,  the  subsystems  are  affected  by  electromagnetic 
energy  from  outside  the  system  as  well  as  from  inside  the  system  which  will 
be  denoted  as  intersysf.em  effects  and  intrasystem  effects  respectively.  Inter- 
system  effects  result  from  natural  as  well  as  man  made  sources.  Examples  of 
natural  sources  are  cosmic  or  extraterrestrial  sources  and  thunderstorms.  The 
man  made  sources  include  desired  as  well  as  undesired  or  extraneous  energy. 

Hie  radiation  from  neon  signs  and  spurious  emissions  from  a  transmitter  are 
examples  of  undesired  man  made  sources,  whereas  the  fundamental  frequency  emis¬ 
sion  of  a  transmitter  is  an  example  of  a  desired  man  made  source.  It  seems 
reasonable  to  assume  that  intrasystem  effects  are  the  more  difficult  to  con¬ 
sider  if  we  assume  that  the  physical  proximity  of  systems  is  much  greater  than 
that  of  the  subsystems  within  a  particular  system.  In  view  of  this  assumption, 
one  may  model  intersystem  e* feet?  <ts  point  sources  or  uniform  field  sources  in 
many  cases.  Additionally,  in  determining  intersystem  coupling  between  two  dis¬ 
tinct  systems,  the  wire-to-wire  coupling  between  systems  seems  intuitively  to 
be  negligible  in  comparison  with  intrasystem  wire-to-wire  coupling.  Also  an- 
tenna-to-antenna  coupling  for  intrasystem  analysis  will  require  near  field  con¬ 
siderations  in  many  cases  whereas  far  field  considerations  are  generally  suffi¬ 
cient  for  intersystem  considerations. 

In  the  past,  intrasystem  coupling  of  electromagnetic  energy  generated 
within  the  system  has  been  considered  through  the  adjustment  of  the  general 
subsystem  specifications  by  granting  waivers.  This  adjustment  of  subsystem 
specifications  is  based  upon  engineering  experience  and  results  in  many  cases 
in  overdesign  with  associated  unnecessary  costs  and  in  some  cases  in  system 
incompatibilities  necessitating  costly  retrofit.  These  inaccuracies  should 
not  be  totally  attributed  to  the  project  engineers,  In?  physical  size  and 
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complexity  of  modern  systems  present  s  monumental  problem  of  interaction  which 
in  many  cases  exceeds  human  cnalytical  capability, 

Pearlston  has  pointed  out  many  inconsistencies  in  the  primary  subsystem 

ni 

EMC  specifications;  MIL-STD-461,  462,  6181,  and  326  .  The  specifications 

basically  limit  the  amount  of  undesired  radiated  and  conducted  electromagnetic 
energy  from  equipments  (e.g.  cable  conduction,  box  radiation,  etc.)  and  provide 
upper  bounds  on  the  susceptibility  of  a  device  to  radiated  and  conducted  elect¬ 
romagnetic  energy  which  enters  a"  various  "ports"  of  che  device  (e.g.  wires, 
antennas,  etc.). 

The  obvious  problem,  however,  is  not  the  subsystem  specifications  but 
that  the  totality  of  the  coupling  paths  is  not  presently  considered  in  suffi¬ 
cient  analytical  detail.  For  example,  one  could  model  the  transfer  paths  (e.g., 
wire-to-wire,  v  >4.w-wire,  antenna-to-wire,  etc.)  with  mathematical  transfer 
functions  which  accurately  reflect  the  physical  location  and  proximity  of  the 
subsystems  within  the  system  and  the  types  of  coupling  paths  (e.g.,  conducted, 
radiated,  wire-to-wire,  etc.).  The  inputs  to  these  mathematical  transfer 
functions  could  be  the  limits  in  the  above  applicable  MIL-STD's  for  the  gen¬ 
erated  extraneous  (undesired)  energy  (we  are  assuming  as  a  worst  case  that  the 
devices  actually  emit  signals  at  the  levels  of  the  limits)  and  design  deter¬ 
mined  levels  for  desired  energy.  To  consider  the  total  electromagnetic  energy 
within  a  system,  one  must  include  all  desired  sources  of  electromagnetic  energy 
(e.g.,  transmitter  fundamental)  as  well  as  all  undesired  sources  (e.g.,  trans¬ 
mitter  spurious  frequencies).  The  total  effect  of  these  sources  on  each  sus¬ 
ceptible  device  or  receptor  could  be  determined  through  an  in-depth  quantita¬ 
tive  analysis  based  on  the  mathematical  description  of  the  various  coupling 
paths.  The  total  electromagnetic  effect  at  the  receptor  port  could  be  com¬ 
pared  with  the  susceptibility  limit  for  the  appropriate  MIL-STD  to  determine 
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if  system  compatibility  will  be  achieved  by  subsystem  compliance  to  thes% 
specifications.  The  point  to  be  made  here  is  that  *'*  all  sources  of  electro¬ 
magnetic  energy  (desired  and  undesired)  are  considered  in  determining  the 
electromagnetic  energy  incident  upon  a  subsystem  or  device  and  if  all  sources 
of  extraneous  energy  conform  tn  the  above  applicable  MIL-STD's  then  the  system 
will  be  guaranteed  compatible  (within  the  confidence  of  the  models)  if  the 
total  energy  incident  upon  each  subsystem  (determined  analytically  above)  docs 
not  "violate"  the  applicable  susceptibility  specification  in  the  case  of  unde¬ 
sired  receptors  or  a  design  determined  susceptibility  limit  for  desirec  recep¬ 
tors.  This  of  course  assumes  that  all  devices  ,md  subsystem?  will  be  properly 
tested  to  insure  their  compliance  with  the  appropriate  M1L-STD' i.  It  should 
be  clear  that  if  on  in  depth  quantitative  analysis  such  as  the  above  is  not  per 
formed,  then  little  confidence  can  be  placed  in  system  compatibility  from  com¬ 
pliance  to  the  MIL-STD's. 

The  above  approach  has  difficulties  however.  In  order  to  insure  compati¬ 
bility  one  must  consider  all  sources  of  electromagnetic  energy,  undesired  and 
desired,  intrasystem  and  intersystem.  This  will  be  impossible  to  implement 
even  on  the  largest  computers  for  even  moderate  size  systems.  One  then  must 
make  assumptions  (hopefully  reasonable  ones)  to  made  the  computation  tractable. 
The  mathematical  modeling  of  the  generators,  coupling  paths  and  receptors  can 
be  simple  or  difficult  depending  upon  the  effort  one  wishes  to  place  in  this 
area.  The  modeling  of  the  receptors  or  processors  of  electromagnetic  energy 
seems  to  be  the  more  difficult  of  the  three.  The  processing  of  electromag¬ 
netic  energy  by  these  receptors  is  in  many  cases  quite  subjective  and  does  not 
admit  a  simple  description.  For  example  the  level  of  noise  which  produces 
interference  in  a  communications  receiver  is  highly  dependent  upon  the  listener 

(i.e.,  a  novice  listener  or  a  "seasoned  veteran"). 

7 


A  common  (and  quite  reasonable)  assumption  in  the  determination  of  the 
coupling  descriptions  are  the  assumptions  of  linearity  and  time- invariance  of 
the  coupling  mediums.  By  linearity  of  the  coupling  medium  we  mean  the  following. 
Consider  two  sources  of  electromagnetic  energy  and  and  one  receptor,  R^. 
Basically,  linearity  of  the  coupling  paths  means  that  if  produces  effect 
E ^  incident  upon  with  "turned  off"  and  G^  produces  effect  Ej  incident 
upon  with  G^  "turned  off"  then  G^  and  G^  operating  simultaneously  will  pro¬ 
duce  an  effect  +  E2  incident  upon  R^.  Clearly  this  is  not  a  rigorous  defi¬ 
nition  of  linearity  but  should  serve  our  purposes  here. 

Basically  the  assumption  of  time-ir».ariance  of  the  coupling  paths  means 
that  the  mathematical  parameters  sufficiently  describing  the  coupling  paths  do 
not  change  with  time.  These  assumptions  of  linearity  and  time -invariance  of 
the  coupling  paths  merely  simplify  the  mathematical  analysis  and  have  been 
implicitly  assumed  almost  universally  (and  with  gord  justification)  In  the 
field  of  EMC.  With  these  assumptions  we  ma”  determine  frequency  domain  trans¬ 
fer  functions  for  coupling  paths  since  a  frequency  domain  description  of  a 
linear,  time-atvarlant  system  is  equivalent  to  a  time  domain  description. 

The  inputs  to  the  transfer  functions  can  be  frequency  domain  description 
of  the  generator  Outputs.  For  example,  the  time  domain  output  of  a  generator 
can  be  represented  equivalently  in  the  frequency  domain  through  the  Fourier 
Transform.  The  result  of  transferring  the  generator  outputs  through  the  trans¬ 
fer  functions  of  the  coupling  paths  yield  e  frequency  domain  description  of 
the  effects  of  the  generators  on  the  receptor  which  if  desired  may  be  converted 
to  a  time  domain  description  through  the  inverse  Fourier  Transform.  At  this 
point  the  modeling  difficulties  increase  greatly.  What  is  required  is  a  mathe¬ 
matical  description  of  the  effect  of  this  incident  energy  upon  the  receptor. 
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This  is  clesrly  very  subjective  and  difficult  to  provide  correctly  in  most 
caces. 

For  example,  one  aiay  wish  to  determine  a  frequency  dependent  susceptibil¬ 
ity  threshold  for  the  receptor.  A  common  way  of  doing  with  is  to  apply  sinus¬ 
oidal  input  signals  A  sinmt  to  the  input  port  and  at  each  frequency  f  * 
increase  A  until  a  "malfunction"  occurs.  Then  if  the  total  incident  effect  due 
to  all  the  generators  at  the  receptor  port  exceeds  this  threshold  at  a  particu¬ 
lar  frequency  then  surely  one  can  say  that  the  device  is  "susceptible."  However, 
if  the  incident  effect  (described  in  the  frequency  domain)  at  a  particular  fre¬ 
quency  does  not  exceed  the  threshold  at  any  frequency  then  can  one  say  that  the 
device  is  not  susceptible?  It  is  doubtful  that  one  can  be  assured  of  this  for 
the  following  reason.  If  the  receptor  processes  the  incident  "effects"  linearly 
then  we  may  represent  this  as  a  linear  system.  If  an  input  u(t)  is  applied  to 
a  linear,  time-invariant  system  with  impulse  response  h(t),  then  the  output, 
y(t)s  will  be  given  by  the  convolution 
00 

y(t)  -  J  h(t-T)  u(t)  dT  (1-1) 

-00 

If  the  Fourier  Transform  of  u(t),  U (w) ,  is  obtained,  then  the  output  "spectrum" 
will  bo  given  by 

Y (»)  «  H(«)  U(»)  (1-2) 

where  H(<*»)  is  the  Fourier  Transform  of  h(t).  The  inverse  Fourier  Transform  of 
Y(m)  yields  y(t)  as 

09 

y(t)  m  2ff  J  H(»)  U(«)  d*  (1-3) 

and  the  output  is  a  summation  of  the  effects  of  each  component  of  U(o).  There¬ 
fore  the  susceptibility  level  determined  by  applying  sinusoids  to  the  device  in 
the  above  manner  may  not  be  the  rele\ant  indicator  of  susceptibility  due  to  this 
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summation  of  effects  even  for  receptors  which  process  the  incident  effects 
linearly.  This  ides  of  *  frequency  dependent  susceptibility  threshold  has  been 
ii  constant  use  in  the  field  of  EMC  and  clearly  needs  sosie  justification.  Some 
general  atteapt  is  aside  in  the  MIL-STD's  to  account  for  this  effect  through  the 
use  of  a  broadband  specification.  However,  it  seeois  that  the  reason  for  insti¬ 
tuting  this  type  of  specification  was  the  realization  that  measureoient  equip¬ 
ments  do  not  have  infinitesloially  small  bandwidth. 

To  simplify  the  analysis,  one  might  also  make  the  assumption  of  "unilateral 
effects."  By  this  we  mean  that  the  receptors  do  not  "load"  the  generators  and 
one  may  simply  determine  a  "transfer  function"  between  each  generator  and  each 
receptor  whose  output  is  dependent  on  the  input  yet  the  input  is  not  affected 
by  the  output.  Also  one  might  assuise  that  classes  of  generators  do  not  inter¬ 
act.  For  example,  we  might  assume  that  the  output  of  a  transmit'  ..  is  not 
affected  by  the  wire  voltages.  If  these  assumptions  are  made,  then  the  compu¬ 
tational  difficulties  are  reduced  substantially  since  one  my  determine  through 
Independent  analyses  the  effects  at  a  receptor  due  to  each  source  individually. 
Then  however,  all  these  effects  must  be  added  to  yield  the  total  effect  cl  all 
sources  ^n  the  receptor.  This  last  point  is  very  Important  if  one  wi'.nes  to 
discuss  compatibility.  Considering  partial  eifects  of  sources  is  jot  a  valid 
way  of  determining  compatibility. 

Another  difficult  area  in  this  type  of  modeling  effort  would  be  the  deter¬ 
mination  of  intersystem  effects.  However,  as  pointed  out  above,  these  effects 
could  be  modeled  on  a  worst  case  basis  much  more  easily  than  intrasystem  effects. 

The  main  deficiency  in  the  present  approach  system  EMC  is  that  the  energy 
coupled  between  subsystems  is  rot  presently  corsidusd  in  sufficient  analytical 
detail.  Also,  intersystem  effects  are  not  c ur.^idered  except  through  adjustment 
of  the  limits  when  intersystem  sources  ar>  known  (e.g.,  when  a  high  power  radar 


10 


is  in  close  proximity  to  the  system).  We  should  also  point  out  that  the 
MIL-STD's  generally  exclude  functional  signals  from  consideration  (e.g.,  the 
60  HZ  power  frequency).  However,  even  if  compliance  to  the  above  MIL-STD's 
was  shown  to  insure  system  compatibility  through  an  indepth  analysis  as  sug¬ 
gested  above,  there  would  be  no  optimization  in  the  design  of  interference  sup- 
pr  ‘ssion  since  the  MIL-STD's  are  formulated  apriori,  e.g.,  without  consideration 
of  the  needs  of  the  specific  system. 

This  is  not  to  say  that  there  is  no  need  for  general  type  subsystems  speci¬ 
fications  as  are  presently  In  existence.  There  are  many  cases  in  which  the 
final  system  configuration  is  either  unknown  or  not  fixed.  Also  the  system 
environment  may  not  be  known  or  fixed  as  with  an  aircraft.  However,  more  mean¬ 
ingful  general  type  specifications  could  be  developed  for  classes  of  general 
systems  through  the  long  term  quantitative  experience  gained  in  analyzing  par¬ 
ticular  types  of  these  systems  through  an  indepth  quantitative  analysis  based 
on  the  evaluation  of  the  various  coupling  paths  within  each  system. 

A  comprehensive  survey  of  the  problems  and  a  proposed  solution  in  the 
determination  of  Intrasystem  compatibility  can  be  found  in  [2], 

In  conclusion,  we  see  that  the  application  of  subsystem  limits  in  their 
present  useage,  e.g.,  MIL-STD-461,  is  adequate  in  many  cases  for  the  follow¬ 
ing  reasons. 

1)  Simply  placing  limits  on  subsystems  will  not  insure  compatibility 
of  these  subsystems  unless  the  various  coupling  paths  are  con¬ 
sidered  in  detail,  analytically, 

2)  Subsystem  compliance  to  the  specifications  generally  wil)  not 
result  in  optimum  system  design  from  the  standpoint  of  required 
interference  suppression,  minimum  system  weight  and  complexity, 


]  1 


etc . 


3)  The  inpact  of  granting  waivers  on  the  overall  system  compati¬ 
bility  cannot  be  easily  or  accurately  assessed. 

4)  The  inpact  of  future  design  changes,  i.e.,  system  modifications 
in  terns  of  subsystem  additions  or  deletions,  cannot  be  easily 
or  accurately  assessed. 

5)  The  impact  of  design  trade-offs  on  system  compatibility  cannot 
be  easily  or  accurately  assessed. 

6)  The  expense  of  system  tests  generally  exceeds  that  which  is  re¬ 
quired  since  critical  areas  cannot  be  determined  which  would 
allow  focusing  of  the  test  effort. 

7)  Subsystem/equipment  teat  data  cannot  be  related  to  system  per¬ 
formance. 

The  quantitative  analysis  suggested  above  would  provide  a  way  of  determin¬ 
ing  more  precisely  the  cost  and  mission  effectiveness  of  a  system. 
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II.  TECHNICAL  INTRODUCTION 


Here  we  will  discuss  the  concepts  o£  mathematically  modeling  the  transfer 
or  coupling  paths.  We  will  assume  that  the  transfer  paths  consist  of  radiation 
and  direct  conduction  (for  wires  only)  between  wires,  antennas  and  boxes.  We 
then  are  faced  with  modeling  the  coupling  modes  of  wire-to-wire,  wire-to-antenna 
wire-to-box,  antenna-to-wire,  antenna-to-antenna,  antenna-to-box,  box-to-wire, 
box-to-antenna ,  and  box-to-box. 

2.1  Wire-to-wire 

We  will  begin  the  discussion  by  reviewing  basic  transmission  line  theory. 
These  results  and  derivations  are  given  in  many  texts ^  Consider  the 

parallel  pair  of  wires  (e.g.,  a  transmission  line)  of  length  X  and  spacing  be¬ 
tween  centers  D  shown  in  Figure  1.  We  can  model  the  transmission  line  as  in 
Figure  2  where  r^,  r^  are  the  self  inductance  and  resistance  per  unit 
length  for  wires  1  and  2  respectively,  c,g  are  the  capacitance  and  conductance 
per  unit  length  betw'.en  the  wires  and  m  is  the  mutual  inductance  between  the 
wires  per  unit  length. 

The  voltage  and  current  along  the  line  are  functions  of  two  variables  x 
and  t  where  x  represents  distance  along  the  line  and  t  represents  time,  i.e., 
v(x,t),  i(x,t).  Therefore  we  may  write  using  elementary  lumped  circuit  analy¬ 
sis  technique*  for  «  small,  infinitesimal  length  of  line,  Ax 

v(x  +  Ax,  t)  -  v(x,t)  -  -f&  d  ^X|t^  -  rAx  i(x,t) 

(2-1) 

i(x  +  Ax,  t)  -  i(x,t)  »  -cAx  — -  g  Ax  v(x  +  Ax,  t) 
where  the  current  in  line  1,  i^(x,t),  equals  the  current  in  line  2,  ^(x.t), 
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1 5 


FIGURE 


i'.e.,  i^x.t)  -  -i2(x,t>  “  <*nd  1  "  ^*1  +  h.  +  ^  and  r  “  *rl  +  r2^* 

Dividing  (2-1)  by  &x  and  taking  tha  Halt  as  Ax  *  0  yields  the  standard 
partial  differential  equations  describing  the  transaission  line 

(2-2) 

SMiiil  »  -c  .  g  v(,,t) 

0X  0t 

Mow  if  we  assume  all  variables  on  the  line  are  sinusoidal  of  frequency  f  ■  ®/2rr 
then  we  may  write  without  loss  of  generality 
v(x,t)  -  V(x)eJ<Bt 

imt  (2'3) 

i(x,t)  -  l(x)eJW: 

where  V(x)  and  l(x)  are  functions  only  of  x  (this  amounts  to  the  usual  separa¬ 
tion  of  variables  swthod  of  solving  partial  differential  equations). 
Substituting  (2-5)  into  (2-2)  yields 


(^)  'J“  -  -J»‘  I«'J“  - «•  w1* 

(^)  «JW  -  -]•  v(,).i“  -  8  v(x,.JW 

1(St 

By  cancelling  e  we  have 

dV(x)  ■  -  (1»jl  +  r)  I(x)  ■  -  t  X(x) 
dx 

(2-5) 

dl(x)  -  -  ()<*c  +  g)  V(x)  -  -  y  V(x) 

dx 

a  set  of  ordinary  differential  equations.  Differentiating  (2-5)  yields 
d2V(x)  -  -  z  “  *y  V(x) 

dl2  *  (2-6, 

and  similarly 

-  zy  l(x)  (2-7) 

dx2 
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These  equations  are  solved  from  (D  ■  —  operator). 

(D2  -  zy)  V(x)  ■  0 

There  are  two  roots  of  che  characteristic  equation  (D^  -  zy) 
the  solution  is 

V(x)  -  Ae**  +  Be-** 

I(x)  -  Ce**  +  De”^* 
where  y  «  J~ty  . 

Differentiating  the  second  equation  of  (2-9)  yields 


(2-8) 

0,  ±  J zy~  »  and 


-  YCe7*  - 


dx 


YDe 


—y  V(x) 


(2-9) 


(2-10) 


(2-11) 


from  (2-5).  Therefore  the  equations  of  the  line  become 

V(x)  -  Ae**  +  Be"** 

I(x)  — Aje**  i  B_e"y* 

*0  *0 

where  z^  ■  J z/y  is  the  characteristic  impedance  of  the  line. 

Now  to  solve  (2-11)  we  must  apply  boundary  conditions  to  obtain  A  and  B. 
Suppose  the  terminations  for  the  line  are  as  shown  in  Figure  3  where  8  and  8 


are  complex  impedances  and  v  (t)  ■  V  e 

8  M 


jwt 


is  a  complex  sinusoid.  Then  we  have 


from  Figure  3  V(0)  *  V  -  1(0)2  ,  V(I)  »  1(1)2  .  Therefore  we  may  solve  the 

dO  L 


following  equations  for  A  ard  B 
V(0)  «  A  +  B  “  Vc  -  1(0)2 

w»  b 

V(X)  «  AeVX  +  Be”*1  -  1(1)2. 


1(0)  «  A_  -  J_ 


(2-12) 


ZQ  S0 


1(X)  *  A_  yX  B__  -yX 
”* D  8  *  -ir  © 


We  then  may  obtain  v{x,t)  =  V(x)eJU,t  and  i(x,t)  *>  I(x)e^Wt  for  any  point  x  on 


th*  lir.e  and  any  sinusoidal  voltage  source  v  (t)  *  V  e 

8  b 


jwt 
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Zs  1(0)  l(J) 


FIGURE  3 
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We  could  also  produce  a  set  of  equations  relating  1(0),  V(0),  I(X)  and 
V(X),  the  terninal  variables  of  the  line,  in  two  port  parameter  form  as 

V(0)  -  V(X)  cosh  yX  +  I(X)za  sinh  yX 

0  (2-13) 

1(0)  ■  I(X)  cosh  yX  +  V(X)  sinh  vi 

20 

Now  the  difficult  part  is  computing  the  capacitance,  inductance,  resistance 
and  conductance  per  unit  length  for  the  lines.  To  illustrate  typical  calcula¬ 
tions  for  these  quantities  we  will  consider  a  pair  of  infinitely  long  conduct  >cs 
as  in  Figure  4.  The  wires  have  radii  r  and  r  .  Let  us  assume  that  wire  1 
and  wire  2  carry  q^  and  q^  coulombs  of  charge  per  unit  length  distributed  uni¬ 
formly  over  the  surface.  Consider  a  single  wire  in  free  space  with  q  coulombs 
per  unit  length  as  in  Figure  5.  The  electric  field  due  to  this  charge  distri¬ 
bution  will  extend  radially  from  the  conductor  since  the  charge  is  assumed  to  be 
uniformly  distributed  over  the  surface  of  the  wire.  From  Gauss's  Law  (which 
also  is  one  of  Maxwell's  equations),  "The  surface  Integral  of  the  normal  com- 

»i»V 

pcnent  of  the  electric  flux  density,  D,  0'’zz  any  closed  surface  equals  the 
charge  enclosed,  Q,"  i.e., 

JJ  D  •  da  ■  Q  (2-14) 

where  da  is  a  unit  surface  area  normal  to  D.  If  we  assume  that  the  medium  is 
homogeneous  and  isotropic  then  D  =  cE  where  •  is  the  permittivity  of  the  medium, 
then  (2-14)  becomes 

J ‘J  I  •  A  -  «/*  (2-15) 

Therefore  using  Gauss's  Law  where  the  area  is  a  "pill  box"  with  the  wire  pass¬ 
ing  through  its  center,  of  R  radium  and  length  1  meter, 

Q  *  q  x  (1  meter)  ■  f f  "D  •  da  »  «E  2rrR  (2-16) 

Since  the  surface  area  is  2ttR  x  1.  Therefore,  the  magnitude  of  E,  |E|,  is 
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given  by 

1*1  ‘  2k  <2'17> 

and  ia  directed  radially  from  the  wire. 

Now  looking  "end-on"  on  the  wire  we  will  calculate  the  potential  differ¬ 
ence  between  two  points  P^  and  as  in  Figure  6.  In  calculating  the  potential 
difference  between  and  P^  we  use  the  definition  of  potential  difference 


'12  J 


n 

r  *■ 


d  t 


(2-18) 


where  df  ia  the  unit  path  length.  Along  the  integral  is  zero  since  this  is 
an  equipotentlal  surface  (due  to  the  fact  that  E  •  df  «*  |eJ  |df|  cos  8  and  9 
along  this  path  is  90°).  Along  path  we  simply  have 


'12  "  i 


q  d  £ 
2n«f 


(2-19) 


—Jr  £n  ~  volts 
2rr*  Dx 


The  capacitance  between  two  bodies  is  defined  as  the  ratio  of  charge  on 
the  two  bodies  to  the  potential  difference  between  the  two  bodies,  i.e., 

C  -  9/v  (2-20) 

The  capacitance  per  unit  length  of  the  two  wire  line  is  therefore 


c  -  q/V_, 


(2-21) 


where  q  is  the  sum  of  the  charge  per  unit  length  for  each  line  and  v^  *s  the 
potential  difference  between  the  conductors. 

Returning  to  the  original  problem  of.  the  determining  the  capccitance  be¬ 
tween  the  two  vires,  let  us  assume  the  "medium”  is  linear  and  use  superposition. 
Assume  qj  *  0.  Therefore  wire  2  being  a  conductor  is  an  equipotential  surface 
and  the  potential  difference  due  to  is  from  our  previous  derivation 


!l_  JL 

rml 


2"f 
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(2-22) 


equipotential 
surfaces  — 


FIGURE  6 
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where  we  use  the  path  of  integration  shown  in  Figure  7,  and  and  are 
equipotential  surfaces. 

Now  we  find  due  to  q^  with  »  0  as 


-  q?  £n  D  ^  q2  In  r»2 
2n«  r##2  2n«  D 

Therefore  the  total  potential  difference  is 

{*£  +  *■¥} 

where  we  assume  q^^  *  -q  *  q.  Rewriting  (2-24) 

,2 


12  2tt« 
and  therefore 


»(  —  ) 

'  r#»lr®2  ' 


(2-23) 


(2-24) 


(2-25) 


2rr« 


'w  » (r®F) 


(2-26) 


Now  we  compute  the  inductance  per  unit  length,  l.  Consider  a  single  wire 
in  a  medium  of  permeability  p.  Vo  assume  a  current  1  is  uniformly  distributed 
throughout  the  wire.  Inductance  l  (per  unit  length)  is  defined  as  V I  where  (r 
is  the  total  amount  of  flux  (per  unit  length)  linking  (enclosing  in  a  complete 
path)  the  current  1.  So,  t  is  what  we  are  seeking.  First,  we  determine  tint» 
the  flux  internal  to  the  wire,  and  then  we  determine  i  .  the  flux  external  to 
the  wire. 

Consider  Ampere's  Law 

„  _i  — * 

!  H  •  dA  “  I  (current  enclosed  by  path)  (2-27) 

closed 
path 


how  the  field  in  the  wire  is  totally  tangential  (assnmed  infinitely 
long  wire  and  uniform  current  distribution),  therefore 

J  H  •  di  -  2nr  HT  -  Ir  (2-28) 

where  is  the  tangential  magnetic  field  and  1^  is  the  portion  of  1  enclosed 
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in  radius  r  as  shown  in  Figure  8.  If  we  assume  Che  current  is  uniformly  dis¬ 
tributed  in  the  wire,  then  the  portion  of  I  enclosed  in  radius  r  is  given  by 


I 


r 


(nr2) 


I 


(2-29) 


Therefore 

2 

2nr  H_  *>  ~~T-  I 

"T  "5 

or 

Vsf  1 

and  the  magnetic  flux  density  (per  unit  length)  is 


(2-30) 


(2-31) 


bt  ■  "A  ■  -^5“  1  (2‘32) 

0)1 

where  ^  is  the  permeability  of  the  wire.  In  a  tubular  element  of  thickness  dr 
(0 

ind  unit  length  along  the  line,  the  total  magnetic  flux  (per  unit  length)  is 
the  product  of  B^,  and  the  cross  sectional  area  of  the  element  normal  to  the 
flux  lines  or  $  ■  J  J  B  •  5a  (da  is  the  unit  surface  area)  or 


l*u,r 

d*  -  I  dr 

2nr2, 

(1)1 


(2-33) 


The  flux  linkage  df  per  unit  length  due  tc  the  flux  in  the  tubuAar  element 
is  given  by  the  product  of  the  flux  per  unit  length  and  the  fraction  of  the 


current  linked.  Thus 
2 

wr 

<»♦ 


Ttr 

TTr, 


0)1 


d$  »  '  dr 

2nr> 


(2-34) 


(1)1 


Integrating  from  the  center  of  the  conductor  to  its  outside  edge  to  find  the 


total  flux  linkages  inside  the  conductor, 
^int 


rr«l  Sair?  . 
Jc  2rr4  ar 
Oil 


(2-35) 


8n 
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Now  we  vast  compute  the  flux  linkages  fro*  flux  external  to  the  conductor. 


The  H  field  is  (tangential  to  the  conductor)  and 

*1  “  2nr 

B  «*£- 
T  2nr 


(2-36) 


where  the  permeability  of  the  medium  surrounding  the  wire  is  p.  The  flux  in 
the  tubular  element  of  thickness  dr  i^  (♦  *  J  T  b  •  da)  (aee  Figure  9) 

(2-37) 


■  IS?  dr 


The  total  flux  linkages  between  and  P2  then  is  given  by 

Do 


Text 


*  2  1*1  ,  ,  ul  °2 

JDi 


(2-38) 


Now  consider  the  two  wire  line  in  Figure  10  with  associateo  currents  Ij 
and  directed  into  the  paper.  Let  us  determine  the  flux  linkages  for  wire  1 
with  I2  *  0  between  the  wire  and  point  X.  From  equation  (2-35)  and  (2-38)  we 
have 


.  .  *  ah 

*n  ir  +  —  to  — 


2TT 


ml 


Now 


for  I2  ♦  0  and  1^  ■  0  we 


have 


.  _  *Z2  °2X 

*12  2TT  *  ' 

Therefore  the  total  flux  up  to  point  X  linking  1^  is  given  by 


♦l  "  *11  +  *12 
*1 


Vl 


-  D1X  .  t*I2  .  P2X 

1T  +  2T  1,1  7T  +  W 

ml 


Equation  (2-41)  may  be  rewritten  as 

1  1*I2  1  u  T  1 

*l“'t”+2Tin^+2rin5  +  2?[Il  *D1X  +  I2  D2x] 

Since  "  “I|  (we  may  assume  a  closed  circuit)  we  may  write 
Uji  Ul,  ,  Ul0  ,  „  r  D. 

*1 


8rr 


^  1  ,  j  ,  m.  r,  ix] 

+  2Tln7^  +  7rln5  +  ^L1iln^J 


(2-39) 


(2-40) 


(2-41) 


2X' 


(2-42) 


(2-41) 
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Allowing  point  X  to  go  to  infinity  determine  the  total  flux  linking  conductor 

D1X 

1  due  to  1^  and  1^  aunt  that  ^2X  4  *  £**e  tern  in  (2-43)  vaniahea. 


Thero5ore  we  nay  interpret  the  self  inductance  of  wire  1  as 
u  1 

and  the  uutual  inductance  as 


u  .  1 

”21  -  2?  >“  D 

Similarly  for  the  second  wire  we  have 

#  *  !itt  +  Ji—  |n  —— 

H  Sir  2n  rot  , 
w2 


*12  ‘  ^  D 


(2-44) 


'2-45) 


(2-46) 


Define  a  *  “  ®2l* 

We  may  also  derive  the  per  unit  length  capacitance  to  ground  for  a  wire 
suspended  above  a  ground  plane  as  shown  in  Figure  11,  To  solve  this  problem 
we  make  use  of  the  method  of  images.  Essentially  this  method  states  that  we 
may  effectively  replace  the  ground  plane  with  the  image  of  the  conductor  as  in 


Figure  12.  Fiom  equation  (2-26)  we  obtair 


»  ‘*TTt 

'g  *  fn  4hz 

7T 


TTS _ 

jtn  2h 
r . 


(2-47) 


tt>  0) 

Similarly  one  may  determine  the  self  inductance  of  a  wire  above  e  ground 
plane  through  the  method  of  images.  Also  the  mutual  inductance  and  capacitance 
between  two  wires  above  a  ground  plane  may  be  determined  through  the  method  of 
images , 

The  resistance  per  unit  length,  r,  can  be  written  as 


r  »  r .  +  r  (f) 

d.c.  a.c. 


<2-4«) 


where  r .  is  the  bulk  resistance  of  the  wire  aad  r  (f)  is  a  frequency 
u i c  #  a.c. 

dependent  resistance  reflecting  the  fact  tnat  for  higher  frequencies,  the 
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current  tends  to  concentrate  more  closely  to  the  outer  surface  of  the  conductor, 

i.e.,  skin  effect.  The  above  per  unit  length  line  parameters  are  derived  in 
r  3-8"1 

many  texts.  Thus  we  nay  draw  the  equivalent  circuit  for  an  infinitesimal 

length  of  line  AX  as  in  Figure  13. 

The  derivation  and  solution  of  the  partial  differential  equations  of  the 
line  are  finite  straightforward.  The  difficult  portion  of  the  analysis  is  the 
derivation  from  elementary  field  theory  of  the  per  unit  length  parameters  for 
actual  lines  (finite  length,  shielding,  grounding,  etc.).  Also  in  the  above 
derivations  v.  have  used  a  number  of  contradictory  assumptions  to  make  the  com¬ 
putations  tractable.  For  example,  we  assumed  that  the  current  through  the  wire 
was  uniformly  distributed  throughout  the  wire  in  determining  the  inductance  of 
the  wire  per  unit  length.  However,  the  current  for  higher  frequencies  tends  to 
concentrate  more  closely  to  the  outside  surface  of  the  conductor  due  to  skin 
effect.  Also  the  equation  for  the  capacitance  per  unit  length  given  by  (2-26) 
is  only  accurate  for  D  »  r^  and  D  »  r^.  This  is  due  to  the  fact  that  closely 
spaced  wires  produce  a  nonuniform  charge  distribution;  "proximity  effect."  A 
more  accurate  formula  is  given  in  [5]  and  f7\  Also  we  should  point  out  that 
the  above  calculations  for  self  and  mutual  inductance  per  unit  length  and  capaci¬ 
tance  per  unit  length  were  made  assuming  static  conditions!  Maxwell's  equations 

—ik 

for  the  time  varying  field  case  modify  Ampere's  Law  as  ^)Tl  •  <Tf  «*  J  (j  +  Ip)  •  da 

closed  path  .-^closed  surface 

where  the  term  J  is  the  usual  conduction  current  density  and  ^  is  the  displace¬ 
ment  current  density.  For  static  fields,  the  equation  reduces  to  Ampere's  Law. 
For  non-static  conditions,  i.e.,  sinusoidal  excitation  of  the  line  especially 
for  high  frequencies,  one  thould  be  very  careful  to  determine  the  correct  per 
unit  length  parameters  and  not  assume  that  static  parameters  can  be  used  to 
describe  non-static  conditions.  Despite  these  approximations  and  inconsistencies 
msny  very  accurate  results  are  obtainable.  Actually  one  could  determine  the  per 


i,(x,t)  Ax 

«-* — foT*— 


unit  length  parameters  through  measurements  on  representative  lines  using  the 
above  mathematical  model  without  resorting  to  a  strict  field  theory  derivation. 
Ihe  derivation  of  these  quantities  for  actual  lines  is  a  very  difficult  matter 
(if  one  does  it  properly). 

A  common  way  of  approximating  the  distributed  parameter  model  of  a  two  wire 
line  is  the  following.  If  I  is  the  total  length  of  the  line,  then  define 
C  •  c  X  f ,  LA  -  X  l,  Rj  *  r^  x  I,  G  -  g  x  I,  and  M  ■  m  X  X  which  yield  the 
lumped  equivalents  given  by  the  lumped  T  equivalent  in  Figure  14a  and  the  lumped 
tt  (pi)  equivalent  in  Figure  14b.  One  can  achieve  better  approximations  by  cas¬ 
cading  the  above  sections  to  yield  the  lumped  double-T  equivalent  in  Figure  15a 
and  the  lumped  double-n  equivalent  in  Figure  15b.  One  could  continue  cascading 
these  lumped  sections  and  in  the  limit  as  the  number  of  sections  approach  infin¬ 
ity  we  would  have  the  exact  distributed  model.  (Notice  that  the  total  line 
parameters,  L^,  R^,  G,  C,  M,  are  proportioned  among  the  number  of  lumped  elements 
in  each  representation.) 

Let  us  apply  the  lumped  equivalents  to  the  network  of  Figure  3.  Using  the 

single  section  lumped  TT  network  in  Figure  14b  we  have  the  lumped  equivalent  in 

Figure  16.  In  Figure  16  we  have  a  simple  lumped  network  to  solve  for  the  ratio 

VL/ 

ll(uo)  ■  'Vg.  We  may  use  standard  lumped  network  theory  to  solve  this  problem. 

(We  may  write  Node  or  Loop  equations  or  we  may  assume  ■  1  volt  &L  and  work 

backwards  to  obtain  Vc).  Note  that  with  this  method  we  do  not  need  to  apply 

w 

boundary  conditions  to  determine  A  and  B  in  (2-11)  as  in  (2-12).  However,  this 
method  is  an  approximation  to  the  exact  distributive  model. 

the  crucial  problem  here  is  the  application  of  these  concepts  to  large 
networks  of  wires  and  their  coupling  affects.  What  must  be  kept  in  mind  is  the 
trade  off  between  computational  difficulty  and  accuracy  of  model  representation. 

A  fundamental  assumption  that  seems  to  be  necessary  is  that  all  wire  voltages 
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and  current*  are  sinusoidal.  Without  going  into  the  details  of  the  justifica¬ 
tion  we  aay  point  out  that  the  equations  for  wire  coupling  are  linear  and  time- 
invariant.  Maxwell's  equations  are  linear  partial  differential  equations  if 
the  permittivity  and  permeability  of  the  medium  (e  and  |*)  do  not  depend  on  the 
electric  field  intensltv,  E,  and  the  magnetic  field  intensity,  H,  respectively. 
Therefore  a  frequency  domain  description  is  equivalent  to  a  time  domain  descrip 
tion  through  the  Fourier  Transform  since  superposition  (for  linear  circuits) 
holds.  This  simplifies  the  required  mathematics  considerably.  However,  as 
shown  above  we  still  will  be  faced  with  the  solution  of  ordinary  (instead  of 
partial)  differential  equations  (see  equation  (2-5)). 

A  common  and  quite  reasonable  approximation  to  circumvent  the  need  for 
differential  equations  is  the  use  of  lumped  tt  or  lumped  T  equivalent  networks 
for  re'  esentlng  these  lines  and  their  coupling  effects.  As  we  are  aware  of, 
lumped  networks  yield  accurate  results  only  for  frequencies  where  the  actual 
circuit  dimensions  are  "electrically  short,"  i.e,  for  dimensions  much  less 
than  a  wavelength.  These  lumped  equivalents  hav*  been  used  for  the  past  50 
years  in  the  Power  Industry  but  there  the  frequency(s)  of  interest  is  merely 
60  H3  and  the  regions  of  applicability  of  the  lumped  models  can  be  determined 
accurately.  In  our  case,  however,  we  will  be  considering  frequencies  up 
to  and  possibly  above  the  GHZ  range.  Therefore  for  higher  frequencies,  one 
could  use  the  double  TT,  triple  tt,  etc.  to  give  more  accurate  results.  Keep  in 
mind,  however,  that  these  lumped  equivalent  networks  are  approximations  to  the 
exact  distributive  analysis.  To  provide  a  correlation  between  experimental 
results  and  mathematically  predicted  results  one  must  be  careful  to  use  the 
appropriate  model.  For  low  frequencies,  the  lumped  equivalents  may  yield  gcod 
results.  For  high  frequencies,  it  may  be  necessary  to  use  the  exact  distri¬ 
butive  model. 
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An  exact  analysis  of  coupling  for  multi-conductor  lines  can  be  made  in  the 
following  straight  forward  manner  using  an  exact  distributive  model.  Consider 
the  set  of  n  wires  shown  in  Figure  17  all  coupled  over  a  distance  £,  i.e.,  all 
lines  are  parallel  over  the  distan-e  f.  Suppose  each  wire  has  a  self  inductance 
and  resistance  per  unit  length  denoted  by  i ^  and  r^  respectively.  Suppose  wire 
i  has  mutual  inductance  per  unit  length,  m^  with  each  of  the  other  wires  and 
mutual  capacitance  and  conductance  per  unit  length  c  and  g^  between  wire  i 
and  each  of  the  other  wires.  We  then  may  draw  the  equivalent  circuit  for  an 
infinitesimal  length  Ax  as  shown  in  Figure  18.  The  line  voltages  ^j^***)  am* 
line  currents  I^x.t)  for  i,j  =  l,...,n  i±j  are  defined  in  Figure  17. 

Let  us  assume  that  all  line  voltages  and  currents  are  sinusoids,  i.e., 
V^x.t)  *  Vij(x)e^lUt  and  ^(x.t)  -  Ii(x)e^U)t.  Therefore  we  may  write 

V  (x+Ax)  -  V  (x)  =8..  Ax  1.  (x)  +  ...  +  2  fix  I  (x)  +  ...  +  2  Ax  I  (x) 
ij  ij  il  1  ii  i  in  n 

I^x+Ax)  -  l.(x)  =  YU2  Ax  V12(x+Ax)  +  ...  +  Yi(n_1)n  V<n  l)n  'x- Ax) 
Dividing  by  Ax  and  letting  Ax  <4  0  we  have 


dVii(x) 


dlt(x) 


-  ailI]L(x)  +  ...  +ainin(x) 


(2-50) 


Yil2  Vi2(x)  +  •**  +  Yi(n-l)n  V(n-l)n  (x) 


where  the  a  and  Y,  _  are  complex  impedances  which  are  functions  of  <v  where 
ij  1WB 

u)  ■  2nf  and  f  is  the  frequency  of  interest.  These  quantities  are  obtainaole  in 
a  tedious  but  straight  forward  manner  from  Figures  17  and  18. 

Equation  (2-50)  leads  to  a  set  of  first  order  differential  equations  in 
matrix  notation  as 


(n-1)  (n-1) 


(2-51) 
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FIGURE  16 


FIGURE  17 
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where  .0.  is  the  i  x  j  zero  matrix  with 

i  J  J 

Vx) 

,  V23<X> 

»--v  »-  : 

i  ,  <L  t  : 

-  dx- 

— 


Vx> 


(2-52) 


and  the  (n-1)  x  n  matrix  3  and  the  n  X  (n-1)  matrix  Y  are  given  by 


(n-l)l . a 


(n-l)n 


l(n-l)n 


n(n-l)n 


(2-53) 


It  is  well  known  from  the  theory  of  state-variable  equations  that  the  general 
solution  to  (2-51)  is  given  by 


1  <x)  '  !ie  + 


*2n-lx 


V2n-lX 


(2-54) 


1  <X>  ■  il'  +-+i2n-l* 

where  the  (n-1)  x  1  vectors  V4  and  the  n  x  1  vectors  1^  for  i  *  l, 
2n-l  are  arbitrary.  The  (2n-l)  complex  numbers  y^  1  ■  l,...,2n-l  are 
eigenvalues  and  are  determined  from 


(n-1)  (n-1) 


ki  2n-l 


0 

n  ~  n 


(2-55) 


where  ^n-l  *8  (2n-l)  X  (2n-l)  identity  matrix^^  . 

To  evaluate  these  constant  vectors  in  (2-54),  we  apply  boundary  con¬ 
ditions  determined  by  the  termination  networks  and  to  obtain  Vj^O)* 

V« . Vo.(0)'  Vi)«(I)-  'l®1’  V» . V0)>  V”- 
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There  are  various  computerized  methods  of  determining  the  above  eigenvalues 


and  vectors. 


The  above  type  of  solution  for  multiconductor  transmission 


lines  is  a  natural  extension  of  the  two-wire  case  which  we  discussed  previously 

.  .  .  .  .  .  f  [12-16,42-47] 

and  has  been  used  in  various  forms 

If  this  exact  (within  the  confidence  of  the  per  unit  length  line  parameters) 
analysis  is  judged  to  be  too  difficult  or  consumes  too  much  computational  time, 
then  an  approximate  lumped  analysis  can  he  used  vis  the  lumped  n  (pi)  or  lumped 
T  networks  with  various  numbers  of  stages  in  cascade  as  in  Figures  14,  15  and 
16.  Keep  in  mind  that  this  will  be  an  approximate  analysis  especially  for  higher 
frequencies. 

There  are  many  other  approaches  to  cable  coupled  interference  which  use 

various  assumptions.  In  general,  these  approaches  use  lumped  circuit  concepts 

for  coupling  between  wire  pairs  and  give  valid  results  only  for  low  frequencies 

T  20-23] 

(or  equivalently  "electrically  short"  circuit  dimensions) 


2.2  Field-to-wire  and  wire-to-field 

Obviously,  wires  will  be  susceptible  to  electric  and  magnetic  fields  gen¬ 
erated  by  antennas  and  metallic  boxes.  What  is  needed  is  a  model  for  determin¬ 
ing  this  vulnerability  of  wires  and  the  effect  of  wire  radiation  on  antennas  and 
boxes.  We  have  considered  and  modeled  the  effect  of  wire  radiation  (or  coupling) 
on  other  wires.  Since  antennas  and  metallic  boxes  emit  and  are  susceptible  to 
electric  and  magnetic  fields,  it  is  sufficient  to  determine  the  fields  produced 
by  wires  and  the  «'  sceptibi lity  of  wires  to  external  fields. 

2.2.1  Field-to-wi re 

Let  us  fir*t  consider  the  field-to-wire  problem.  Let  us  suppose  that  we 
have  e  pair  of  wires  (a  transmission  line)  with  an  incident  electric  and  mag¬ 
netic  (ield  as  shown  in  Figure  19.  We  will  give,  two  derivations  of  the  equations 
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FIGURE  19 


FIGURE  20 


which  relate  the  incident  fields  to  the  voltages  and  currents  on  the  lines.  The 
first  is  given  in  [24],  This  derivation  takes  into  account  the  distributed 
nature  of  the  transaission  lines  and  is  exact  in  this  sense.  The  second  deriva¬ 


tion  uses  simple  field  relations  and  lumped  circuit  concepts.  Other  allied 
results  appear  in  [25-27].  A  great  deal  of  interest  in  the  field-to-wire  coup¬ 
ling  problem  is  generated  by  the  need  for  determining  BMP  effects  on  electronic 
circuitry. 


First,  consider  Maxwell's  equations  in  point  form 

V  XK-f+f 

V  xt 

V  dt 


[7] 


(2-56) 


D 


V 

V 


e 

0 


Stokes  Theorem  states  that 
n  -5.  a 

J  ^7  A  *  da  •  §  A  *df 

surface 


(2-57) 


l.e.,  integrating  the  curl  of  a  vector  quantity  over  a  surface  is  equivalent  to 
obtaining  the  line  Integral  of  the  vector  quantity  along  the  closed  perimeter 
enclosing  the  surface.  The  direction  of  the  line  integral  is  such  that  the 
direction  of  da  complies  with  the  right  hand  rule.  Integrating  (2-56)  in  accor¬ 
dance  with  (2-57)  produces  Maxwell's  equations  in  integral  form  as 

A 


4  (f  +  3 


da 


$  *  •  «  -I 


dB 

at 


(2-58) 


da 
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(2-58) 


J  f  -  S?- J  ^ 

Closed  Vol 

Surface 
Enclosing 
Volume 

-i 

J  B  •  da  ■  0 

Closed 

Surface 

Enclosing 

Voluate 

where  the  last  two  equations  of  (2*58)  follow  from  (2*56)  and  Gauss's  Divergence 
Theorem,  i.e., 

J  y'  •  A  dv  ■  J  A  •  da  (2*59) 

Volume  Closed 

Surface 
Enclosing 
Volume 

Consider  a  differential  length  of  line  As  as  shown  in  Figure  20  where  the 
per  unit  length  parameters  of  the  line  are  included.  Assume  without  loss  of 
generality  the  incident  fields  are  sinusoidal  and  consist  of  components 


E^x.z.t)  -  Es(*,*)e 
E£(x,z,t)  -  E£(x,z)e 


j(Dt 

Jaat 

jUJt 


(2-60) 


By(x.z.t)  -  By(x,z)e 
Integrating  the  second  equation  of  (2-58)  yields 
D 

J  [B  (x.z-h&z)  -  E  (x,s)3  dx 
0 

z4Az 

-  J  [E£(D,z)  -  E£(0,z)j  dz 
z 

z+tiz  D 

•  -ju>  J  J  By(x,z)  dx  dz 


(2-61) 
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tfc  aay  easily  identify 


11a 

Az  «♦  0 


(z)  ,  dV(z)  lia  1  „D  r  , 

dz  "  Az  0  fiz  J  *  E,(*»*)]  d* 


lls  1 

Az  -*  0  AZ  I  [E*  <D»*>  '  E*  <°»z>  3  d*  -  ai(*) 
Z> 

where  2  ■  Jj  +  2^  and  therefore  have  the  relation 

*  *I(z)  -  Jut  J  B  (x,zj  dx 
0  y 

Siailarly  we  My  obtain 

&{&  +  YV(z)  -  -  Y  Ex(*,z)  dx 
Frca  .heae  equations,  one  My  obtain 


<V«>  + 

o.2 

and  similarly 


u 

BYV(z)  «  8Y  J*  E  (x,z)  dx 


•> 

-  BYl(z)  -  -  Y[E  (D,z)  -  E  (0,z)  ] 


(2-62) 


(’-£3) 


(2-64) 


(2-65) 


(2-66) 


Equations  (2-65)  and  (2-66)  My  be  solved  in  standard  fashion  after  lnte- 
gnting  the  aource  field  (x ,  z)  and  each  solution  will  contain  two  undeterained 

constants  in  the  homogeneous  part  of  the  solution.  These  constants  My  be  evalu¬ 
ated  with  the  boundary  conditions  V(0),  V(I),  1(0)  and  1(1)  once  one  specifies 
the  terminations  for  the  line.  The  result  will  yield  V(z,t)  »  V(z)eJu,t  and 
I(z,t)  -  l(z)eJtBt. 

A  rather  crude  approximation  can  be  obtained  in  the  following  way.  This 
solution  neglects  the  scattering  effect  of  the  transmission  line.  We  know  that 
a  changing  Mgnetic  field  will  induce  potential  difference  in  «  loop  which  is 
proportional  to  the  rate  of  change  of  flux  enclosed  by  the  loop,  i.e., 
d*  1  D 

Vm  "  ■  dt  "  ’j“>  J*  5  By<**z>  d*  dz  (2-6') 

0  0 

««■  uas  the  lumped  n  equivalent  circuit  for  the  transmission  line  and  include 

Vm  in  the  loop  as  shown  in  Figure  21.  In  Figure  21,  the  per  unit  length  parameters 
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2.2.2  Wire-to- field 

Here  we  will  only  discuss  the  concepts.  It  is  well  known  that  a  current 

carrying  wire  will  produce  a  magnetic  field.  7  .is  magnetic  field  can  be  related 

[6-£ 

to  the  current  through  the  3iot- Savart  Law  '.  Also  the  potential  difference 
between  two  wires  of  a  transmission  line  or  between  two  points  along  a  transmission 
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line  will  give  rise  to  a  radiated  electric  field.  Various  approximations  may 
be  used  to  model  the  transmission  line  as  a  wire-type  antenna  to  obtain  the 
resulting  radiated  fields.  The  exact  solutions  for  radiating  wires  are  well 
known  once  one  is  given  the  current  distribution  along  the  wire.  This  piece  of 
informstion  is  usually  the  s»st  difficult  to  obtain  in  the  solution  of  antenna 
problems. 

Suppose  we  are  given  the  source  and  load  terminations  for  a  two-wire  trans¬ 
mission  line  as  in  Figure  3.  An  exact  solution  for  the  current  on  each  wire  as 
a  function  of  distance  along  the  wire  can  be  found  in  a  sti  \ight forward  manner 
as  given  by  (2-11).  Consequently,  we  may  solve  any  multiconductor  transmission 
line  problem  for  the  currents  along  each  wire  and  use  relatively  straightforward 
integral  field  equations  to  derive  the  spatial  electric  and  magnetic  fields.  He 
then  may  superimpose  the  fields  due  to  each  wire  separately  at  each  point  in  space 
to  obtain  the  electromagnetic  field  profile  due  to  these  radiating  wires.  In 
theory  then,  this  problem  is  quite  straightforward  in  solution.  In  practice,  one 
would  probably  approximate  the  solution  at  desired  points  in  space. 

2.3  Antenna  radiation  and  reception 

If  we  determine  the  electromagnetic  fields  produced  by  an  antenna  then  we 
may  determine  the  effect  of  this  source  upon  wires  through  the  analysis  in  2.2. 

At  this  point,  we  may  Invoke  reciprocity  to  determine  the  reception  capabilities 
of  an  antenna  for  incident  fields.  To  simplify  matters,  we  will  discuss  antenna- 
to-antenna  coupling  through  some  methods  which  hevi'  been  used  a  great  deal  in  the 
past. 

The  exact  solution  for  field-to-antenna,  antenna- to-fleld  and  antenna-to- 
antenna  coupling  has  only  been  solved  in  general,  closed  form  for  filamentary 
wire-type  antennas  (loop,  dipole  and  long  wire).  Consequently,  more  general 
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antennas  yield  solutions  aore  easily  on  an  "input-output"  basis. 


For  deteraining  the  coupling  between  two  antennas,  we  could  assume  as  a 
highly  simplified  but  often  used  approxiisatlon  that  both  antennas  are  isotropic 
point  source  radiators  as  thown  in  Figure  23.  With  this  assumption,  the  radiated 
power  is  given  by 

Pt/4*R2  -  P  (w#tta/«2)  (2-71) 

4  2 

where  P  is  the  power  per  unit  area  passing  through  a  sphere  of  area  4nR  and 
the  transmitter  is  producing  watts  of  power.  In  this  case  the  Poyntlng  vec¬ 
tor  becomes 

P  -  E  x  H  (2-72> 


and  since  we  are  assuming  plane  waves  in  free  space  (i.e.,  isotropic  point  source) 

£|  -  *  |H|  (2-73) 

where  |E |  is  the  magnitude  of  E  and  z  is  the  impedance  of  the  transmission  medium 
(377  ohms  for  free  space).  Thus  we  have 


PT/4rtD2 


(2-74) 


(2-75) 


as  the  magnitude  of  the  electric  field  in  the  vicinity  of  the  receiving  point 
source  antenna.  We  could  also  convert  this  electric  field  to  an  equivalent 
voltage,  V,  at  the  base  of  the  receiving  antenna  using  another  highly  simplified 

A 

(but  often  used)  approximation  with  an  "effective  antenna  height"  h  as 

A 

V  -  h  ]E|  (2-76) 

Clearly  this  is  a  very  unsatisfactory  model  for  reasons  almost  too  obvious 
to  mention.  The  model  does  not  consider  the  peculiarities  of  the  particular 
antenna,  such  as  beam  pattern,  etc.  Also  the  "effective  height"  of  the  receiving 
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antenna,  h,  can  be  directly  calculated  only  for  antennas  of  very  siaple  geome¬ 
tries.  In  most  cases  this  quantity  is  very  subjective  and  is  obtainable,  at 
best,  through  direct  measurement.  The  model,  however,  is  the  simplest  quantita¬ 
tive  model  of  antenna  coupling. 

Define  the  effective  area  of  an  antenna  as  the  ratio  of  the  received  power 
of  an  antenna  to  the  power  density  of  the  Incoming  wave,  i.e.. 


A  - 


R 

Pin 


(2-77) 


Suppose  we  have  a  transmitting  antenna  of  effective  area  and  a  receiving 

antenna  of  effective  area  A^.  A  commonly  used  formula  for  determining  the 
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received  power  of  an  antenna  is  the  Friis  transmission  formulaL  *  *  •> 


PR- 


Pt  \  \ 

A* 


(2-78) 


where  X  “  2»/f 

and  f  is  the  frequency  of  transmission.  It  should  be  noted  that 
in  the  above  usage,  the  effective  areas  Ap  and  must  be  given  in  the  direction 
of  transmission.  This  formula  has  been  used  extensively  today  even  though  it  was 
determined  in  1946.  The  reason  for  its  success  seems  to  be  its  simplicity  and 
the  very  general  definition  in  (2-69).  The  burden  is  placed  upon  the  user  to 
determine  precisely  what  the  effective  area  of  each  antenna  is,  no  matter  how 
complicated  ,  geometrically,  his  antenna  is. 

In  the  case  of  coupling  between  antennas  collocated  on  an  aircraft,  the  trend 
has  been  to  use  the  Friis  transmission  formula  and  a  corr«,:tion  factor  to  take 
into  account  the  bending  of  the  waves  as  they  follow  the  curvature  of  the  fuse¬ 
lage^®’  Wing  shading  and  diffraction  effects  are  included  in  [28], 

For  propagation  around  none yl indr leal  and  especially  unsymmetrlcal  bodies,  the 
analysis  could  be  amenable  only  to  experimental  determination. 
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2.4  Box  radiation  and  wotion 


Very  little  has  been  done  in  Modeling  fleld-to-box,  box-to-field  end  box- 
to-box  coupling,  end  it  ie  doubtful  that  any  aignificant  reaulta  will  be  obtained 
in  the  near  future  due  to  the  variability  in  geoawtry  and  location  of  the  boxes 
(which  are  not  designed  for  functional  coupling  purposes  as  with  antennas)  with¬ 
in  a  system.  This  area  of  coupling  determination  seems  to  be  the  sa>st  difficult 
to  model.  This  is  due  to  at  least  two  reasons.  The  shapes  of  boxes  containing 
electronic  equlpsMnt  is  not  designed  for  functional  purposes  and,  consequently, 
meager  standardised  information  is  available.  Furthenaore,  what  is  required 
essentially  is,  first,  the  determination  of  the  coupling  to  and  from  the  box  wall 
and  external  fields  and,  more  importantly,  the  coupling  between  the  box  and  its 
Internal  electronic  circuitry  (  which  is  the  fundamental  area  of  concern). 

Boxes  containing  400  Hz  transformers  have  been  replaced  by  an  equivalent 
magnetic  dipole  (see  the  discussion  of  the  McDonnell  Program).  Other  allied 
results  are  found  in  [33]. 


III.  THE  BOEING  PROGRAM^3^ 

This  program  was  procured  by  the  Air  Force  Avionics  Laboratory  (AFAL)  to 
analyze  the  intrasystem  EMC  of  avionics  systems  and  is  concerned  with  analyzing 
cable-to-cable  and  direct  cable  conducted  coupling.  To  illustrate  the  program 
we  will  use  a  simple  example.  Consider  the  electrical  network  in  Figure  24. 

In  Figure  24,  the  source  S  represents  a  source  of  electromagnetic  energy 
(either  desired  or  undesired)  and  the  receptor  R  represents  a  device  which  may 
be  critically  affected  by  electromagnetic  energy  (either  desired,  undesired  or 
both).  Let  us  assume  that  these  devices  are  connected  by  parallel  wire  lines 
with  ground  returns  separated  by  distance  D  over  a  common  length  I. 

The  assumption  is  made  that  the  source  S,  load  Z^,  load  and  receptor  R 
are  linear  devices  with  respect  to  their  input  terminals.  This  assumption  is 
Implicit  in  the  representation  of  each  of  these  devices  with  a  Norton  (or 
Thevenin)  equivalent  circuit  as  shown  in  Figure  25. 

In  Figure  25,  Z  (f),  Z.(f),  2-(f)  and  Z_(f)  are  frequency  dependent  impe- 
SI/  K 

dances,  e.g.,  Z(f)  =  R(f)  +  j  X(f)  where  R(f)  and  X(f)  are  the  real  and  imagin¬ 
ary  parts  of  Z(f)  respectively  and  are  functions  of  frequency  f.  The  Norton 
equivalent  current  source,  1(f),  for  S  is  a  frequency  dependent  amplitude  func¬ 
tion  which  is  equivalent  to  the  magnir  •»  of  the  Fourier  Transform  of  the  short 
circuit  current  associated  with  S.  The  phase  angle  of  the  Fourier  spectrum 
associated  with  this  source  is  neglected.  There  may  be  serious  questions  in 
some  cases  about  the  assumption  of  linearity  of  the  devices  (which  is  inherent 
in  the  use  of  a  Norton  or  Thevenin  representation).  However,  the  magnitude  of 
the  problem  seemingly  dictates  this  assumption  as  a  necessity. 
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Receptor 

R 


The  above  Impedances  and  the  spectrum  of  the  current  source  are  determined 
through  direct  measurements  on  the  actual  devices.  The  magnitudes  of  2^(f), 

2^(0  and  2R(f)  given  by  |2^(f)  | ,  ^(f)  |  and  |2^(f)J  respectively  are  deter¬ 
mined  and  an  equivalent  two  terminal  RLC  type  network  is  synthesized  whose  input 

A  A  A 

impedance  are  2j(f),  2^(f)  and  2R(f).  The  synthesized  circuits  are  such  that 

A  A  A 

|21(f)|  -  J21(f)|,  |a2(f ) I  •  J22(f)|,  |zRCf)|  "  !aR(f)|  Where  *  denotes  "approxi¬ 
mately  equal  to."  The  reason  for  determining  an  equivalent  RLC  network  instead 
of  using  the  plots  of  input  impedance  magnitude  vs.  frequency  is  that  the  data 
inputs  to  the  program  are  simplified  since  only  the  values  of  the  elements  and 
the  circuit  configuration  of  the  RLC  equivalent  need  be  specified  in  terms  of 
an  easily  obtained  input  Impedance  equation.  However,  one  could  approximate  the 

A 

curve  |2(f)j  with  straight  line  segments  and  determine  the  equation  for  |2(f)| 
for  this  approximation  to  any  desired  accuracy.  The  equivalent  current  source 
of  the  source,  1(f),  is  simply  a  magnitude  versus  frequency  spectrum.  The  sus¬ 
ceptibility  level  of  the  receptor  is  a  magnitude  versus  frequency  spectrum  of 
the  input  voltage  to  the  receptor,  V„(f),  and  is  determined  external  to  the 
analysis. 

The  transmissions  lines  and  their  coupling  effects  are  modeled  using  the 
single  section  lumped  TT  equivalent  as  in  Figure  26.  Ir.  Figure  26,  L  *  fxi, 

C  *  cxf,  *  c^Xl,  M  ■  ^X£,  R(f)  *  r(f)X£  is  a  frequency  dependent  resistance 
function  which  includes  skin  effect  and  £  is  the  common  length  of  the  generator 
and  receptor  transmission  line. 

What  is  required  is  the  computation  of  the  magnitude  of  VR(f)  to  be  com¬ 
pared  with  the  susceptibility  spectrum  of  the  receptor  to  determine  if  an 
incompatibility  exists,  if  the  magnitude  of  the  received  voltage  spectrum, 

VR(f),  «xceeds  the  receptor  susceptibility  spectrum  at  any  frequency  then  an 

incompatibility  exists.  The  truth  of  the  converge  of  this  statement  is  of 
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course  completely  dependent  upon  having  a  relevant  and  meaningful  susceptibility 
spectrum. 

vD(f)  is  computed  using  the  assumption  that  the  loading  effect  of  the 

K 

receptor  circuit  on  the  generator  circuit  is  negligible.  This  assumption  mani¬ 
fests  itself  in  the  following  way.  The  voltages  V  (f)  and  V .(f)  and  the  current 

a  b 

Ij(f)  in  Figure  26  are  computed  from  the  circuit  of  Figure  27.  With  this  infor¬ 
mation,  the  receptor  voltage,  V  (f),  is  computed  from  the  circuit  of  Figure  28 

K 

where  V  (f)  *  iuMI.(f)  and  uu  =  ?nf. 
ra  1 

The  program  is  applied  to  the  computation  of  interference  for  networks  with 
more  than  one  set  of  emitter  wires  and/or  more  than  one  set  of  receptor  wives  in 
the  following  manner.  Let  us  suppose  that  the  entire  cable  coupling  network  con¬ 
sists  of  p  sets  of  generator  (or  emitter)  wires  G^,  .  G^  and  q  sets  of 

receptor  wires  R, ,  R  ,  ...  ,  R  as  illustrated  in  Figure  29. 

1  2  q 

The  system  we  are  considering  is  then  represented  as  in  Figure  30.  Each 
generator-receptor  pair  is  considered  separately  and  modeled  as  in  Figure  26. 

(The  extension  to  the  case  without  the  ground  plane  is  also  considered.)  A  fre¬ 
quency  dependent  transfer  function  lor  each  generator-receptor  pair  is  deter¬ 
mined  as  described  previously  (see  Figure  26,  27  and  28)  given  by  H^(u>)  where 
the  subscript  i  repre«ents  the  generator  and  the  subscript  j  represents  the 
receptor.  The  reader  should  note  that  this  derivat-on  implicitly  assumes  no 
coupling  between  each  generator,  no  coupling  between  each  receptor  8nd  no  coup¬ 


ling  between  any  generator-receptor  pair  other  than  G.  and  R^  (the  pair  being 
considered).  Furthermore,  it  is  assumed  that  receptor  R^  doe*  not  "load"  gen¬ 
erator  G^  in  calculating  the  transfer  function  H  («b).  In  other  words,  we  are 
deleting  all  generators  and  all  receptors  except  G.  and  R  in  determining  the 

1  j 

transler  function  H  (uj)  .  In  determi  '  ”g  H.  ,  (u))  in  th’  *  reduced  network  we 
J  1 

further  p'-sume  that  the  presence  of  R  does  not  affect  the  values  of  the  voltages 
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ocr 


and  current*  In  the  generator  network  G^.  The  transfer  function  H^(w)  then 
provides 


v* 


Rj 

\  ((0) 
Ii(») 


(3-1) 


under  the  above  assumptions.  These  approximations  were  apparently  made  to  reduce 

computation  time.  To  provide  an  exact  solution  (within  the  confidence  in  the 

line  parameters  and  the  use  of  the  lumped  ft  equivalent),  one  would  have  to  solve 

a  horrendous  network  problem  by  including  all  generator  and  receptor  circuits 

with  all  associated  cross-coupling  effects.  Notice  that  the  voltage  across  the 
G, 

generator  load  2  l(f),  is  calculated  automatically  once  one  solves  the  generator 

L 

network.  Therefore,  the  direct  coupled  interference  transfer  function  is  given 

hy  c. 


VR  <«> 

Ij  (uj) 


(3-2) 


At  this  point,  let  us  discuss  certain  important  concepts  to  be  considered 

in  using  this  program  as  well  as  any  other  compatibility  prediction  program. 

Obviously,  in  the  actual  system  deployment,  the  voltage  arross  a  receptor  irnpe- 
Ri 

dance,  2  (uu) ,  will  be  the  sum  of  the  effect*  due  to  all  generator*  of  inter- 
R 

ference  (not  only  wire  type  generators  but  al*o  antenna,  box  and  intersystem 
sources).  One  cannot  simply  excerise  thi*  program  by  determining  the  voltage 
acro«*  a  wire  type  receptor  load  impedance  due  to  ench  wire  type  source  and 
comparing  each  of  the*e  voltage  spectra  to  the  susceptibility  threshold  spec¬ 
trum  oi  the  receptor  to  determine  compatibility  and  conclude  that  if  none  of  the 
received  voltage  spectroms  exceed  the  susceptibility  threshold  spectrum  the 
‘ystemwfll  be  a  compatible  one,  Thi*  ureage  would  represent  a  gross  conceptual 
error  and  In  many  oa*e*  would  yield  an  incorrect  conclusion.  The  receptor  in 
an  actual  w*tem  will  be  affected  to  some  degree  by  all  source*  (wire,  box, 
antenna  and  intersystem  sources. 


58 


As  an  example  of  the  proper  Method  of  analysis,  let  us  assume  that  the 
entire  system  consists  of  tvo  wire-type  sources  and  one  transmitter  with  asso¬ 


ciated  antenna.  Also  let  us  assume  that  there  is  only  one  receptor  and  this  is 

a  wire-type  receptor.  The  "system"  can  be  modelled  as  in  Figure  31.  In  Figure 

31,  the  coupling  between  the  wire  circuits  represents  mutual  inductance  and 

capacitance  coupling  and  the  coupling  between  the  antenna  and  wire  receptor  is 

only  "one-way."  In  this  model  we  are  assuming  only  that  the  output  of  the 

transmitter  is  unaffected  by  the  rest  of  the  system  and  the  wire  generators,  G, 

and  G are  not  affected  by  the  transmitter.  Also  we  are  assuming  that  there  is 

no  box  radiation  or  reception.  Notice  that  we  are  including  the  effects  of 

Intersystem  sources.  These  intersystem  sources  will  affect  the  entire  system 

and  should  be  included  unless  it  is  assumed  chat  their  effects  are  negligible. 

We  could  then  solve  the  reduced  systems  in  Figure  32,  33,  34  and  35. 

In  Figure  32,  we  simply  "turn-off"  the  transmitter  and  both  current 

sources  I_ (u>)  and  I  (u»)  and  determine  the  contribution  to  V  (u>)  due  to  inter- 
X  A 

system  effects  only,  V  (in)..  In  Figure  33,  we  "turn-on"  the  transmitter, 

Rj  1 

"turn-off"  both  current  sources  Ij(u>)  and  I^tu)  and  intersystem  effects  (in 

making  the  calculations)  to  determine  the  contribution  to  V  (u>)  due  to  the 

R1 

transmitter  only,  V  (u>)  .  Similarly,  in  Figures  34  and  35  we  determine  the 

Rj  2 

contributions  to  VR  (u>)  due  to  the  current  sources  Ij(u>)  and  I^iu)  only. 


respectively. 

Notice  that  if  a  wire  generator  is  "turned-off"  we  do  not  delete  the  wires 

from  consideration.  In  other  words,  the  "deactivated"  wire  generators  still 

exist  but  do  not  have  excitations.  For  example,  in  Figure  33,  the  receptor 

voltage  V  (u>)  due  to  the  transmitter  will  be  affected  not  only  by  the  trans- 
R1  2 

mil (or  output  but  also  bv  the  parasitic  nature  of  the  deactivated  generator 


lines. 
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FIGURE  31 


FIGURE  32 


b'J 


FIGURE  35 


Since  superposition  of  source  effects  is  valid  for  linear  coupling  aedia 


(the  lines  are  linear  and 

the  medium  is  linear). 

the  voltage  across  the  recep- 

tor  load*  V  (t»)f  will  be  given  by 

y>  ■  y*>i + 

Y»>2  + 

VR  (*>3 

_X-< 

+  VR  (»>4 
,  »  , 

(3-3) 
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Now  if  one  compares 

the  received 

threshold. 

VB  (»),  to 
R1 

the  susceptibility 

threshold  of  R^f  a  valid  conclusion  as  to  the  compatibility  of  this  "system"  can 
be  made  (if  the  susceptibility  threshold  of  Rj  represents  the  susceptibility  of 
R^  in  a  valid  manner).  This  is  the  correct  method  of  calculating  the  compati¬ 
bility  of  a  system  and  is  exact  within  the  confidence  of  the  generator,  receptor 
and  coupling  models  along  with  the  confidence  in  the  siodel  of  the  intersyst*m 
effects. 

The  Boeing  program  approximates  this  correct  solution  in  that  the  response 

of  the  receptor,  V  (uo),  would  consist  of  the  sum  of  two  received  spectrums. 

R1 

One  of  which  would  be  due  only  to  Gj  with  the  parasitic  nature  of  G^  neglected 
and  no  consideration  of  the  transmitter  or  Intersystem  effects.  The  other 
would  be  due  only  to  with  the  parasitic  nature  of  Gj  neglected  and  no  con¬ 
sideration  of  the  transmitter  or  intersystem  effects.  If  one  had  "sufficient 
reason"  to  believe  that  the  effects  on  R^  due  to  the  transmitter  end  intersystem 
sources  was  completely  dominated  by  the  effects  due  to  G^  and  G^»  then  within 
this  assumption,  one  should  add  the  spectrums  due  to  G1  and  G0  together.  This 
would  approximate  the  effect  of  G^  and  G ^  on  R^.  It  would  still  be  an  approxi¬ 
mation  since  the  Boeing  program,  for  example,  neglects  the  parasitic  effect  of 
the  deactivated  G^  circuit  when  calculating  the  response  of  due  to  G^..  This 

f>3 
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approximation  of  neglecting  paraaitic  effect*  of  other  generator  lines  (as  well 
as  neglecting  parasitic  effects  of  other  receptor  lines)  as  is  used  by  the  Boeing 
program  is  apparently  felt  to  yield  worst  case  results.  However,  this  may  not 
necessarily  be  true  in  all  cases.  For  example,  resonance  effects  produced  by 
including  the  parasitic  nature  of  lines  other  than  those  of  the  generator  recep¬ 
tor  pair  in  the  computation  conceivably  could  produce  receptor  line  responses  of 
greater  magnitude  than  those  calculated  on  a  generator  receptor  pair  basis  as  in 
the  Boeing  program. 

The  Boeing  program  is  an  extension  of  other  work  in  predicting  cable  coup- 

f  33  37—39 1 

ling  (and  in  certain  programs,  box  radiation  and  reception).  ’  J 
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IV.  THE  McDONNELL-DOUGLAL  PROGRAM L 

This  program  was  procured  by  AFAL  for  the  EMC  analysis  of  airborne  com- 
municat ion-electronic  systems  and  is  currently  in  the  final  stages  of  develop¬ 
ment.  Essentially  the  program  determines  through  analytical  models  of  antennas, 
wires  and  boxes  as  generators  and  receptors  with  associated  coupling  effects 
whether  subsystem  compliance  to  the  appropriate  MIL-STD's  will  insure  compati¬ 
bility  of  the  system  (within  the  confidence  in  the  models  of  the  generators, 
receptors  and  their  coupling  effects).  Functional  as  well  as  nonfunctional 
sources  are  considered. 

The  program  is  modular  in  that  each  of  the  above  coupling  effects  are 
calculated  separately.  For  example,  the  wire-to-wire  coupling  is  calculated 
as  a  separate  routine;  f ield-to-wire  coupling  is  calculated  separately,  antenna- 
to-antenna  is  calculated  separately,  and  box-to-box  coupling  is  calculated 
separately. 

The  wire-to-wire  coupling  routine  uses  the  lumped  rr  equivalent  circuit  to 
represent  the  lines  and  their  coupling  effects.  This  routine  differs  from  the 
Boeing  program  in  that  the  entire  wire  generator,  wire  receptor  network  is 
solved  simultaneously.  In  other  words,  the  computation  is  not  performed  on  a 
generator-receptor  pair  basis  and  the  parasitic  nature  of  other  generator  and 
receptor  lines  is  included  in  determining  the  coupling  between  a  generator- 
receptor  pair.  Nodal  equations  are  written  for  the  entire  network  under  con¬ 
sideration  and  are  solved  by  Gaussian  elimination  for  the  network  node  voltages. 
The  per  unit  length  parameters  are  derived  for  various  shielding,  grounding  and 
twisting  configurations.  For  frequencies  where  the  wire  length  is  greater  than 
1/20  of  a  wavelength  (lor  the  frequency  under  consideration)  the  model  of  the 
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transmission  line  and  its. coupling  effects  consists  of  two  tt  sections  in  c«s~ 
cade.  In  addition,  only  intrabundle  coupling  is  considered  due  to  computer 
storage  limitations  with  the  largest  bundle  consisting  of  approximately  60 
wires  with  4  shields,  40  wires  with  20  shields  or  20  wires  all  with  shields. 

The  calculation  is  performed  in  discrete  frequency  steps  and  the  user  can  input 
the  desired  frequency  range.  The  program  has  the  capability  to  accept  load  and 
source  impedances  for  generator  and  receptor  wires  as  RLC  networks,  since  it  is 
essentially  constructed  as  a  general  network  analysis  program.  The  source  and 
Susceptibility  spectra  are  either  the  appropriate  EMC  specifications,  measured 
spectral  data  (e.g. ,  for  known  functional  signals)  or  an  internal  generation  of 
the  magnitude  of  the  Fourier  Transform  for  the  time  domain  output  waveform  of 
the  generator.  The  program  determines  Intermediate  points  between  the  specified 
discrete  points  in  the  spectrums  by  interpolating  log-linearly  between  these 
points.  The  output  of  the  program  consists  of  the  magnitude  and  phase  of  each 
source  as  well  as  the  total  received  voltage  due  to  all  sources. 

The  antenna-to-antenna  coupling  model  presently  used  is  based  on  the  re¬ 
sults  of  [281.  This  model  is  designed  to  predict  coupling  between  antennas 
collocated  on  an  aircraft.  A  transmitter  antenna-receiver  antenna  pair  is 
selected  and  checked  for  frequency  coincidence.  If  frequency  coincidence  exists, 
then  an  EMI  margin  is  determined  at  the  worst  case  coincident  frequency,  f^, 
through 

EM(f  )  -  PR(f  )  -  SR(f  )  (4-1) 

U)  (C 

where  PR  is  the  power  received  from  the  transmitter  at  the  worst  case  frequency 
and  SR  is  the  sensitivity  threshold  of  the  receiver  at  the  worst  case  frequency. 
The  power  received  from  the  transmitter,  PR,  is  determined  from 

PR<V  -  w +  vv +  w  •  cvv  (4.2) 

-  CW  ■  SF<V  +  20  ">*10  i  &  ] 
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where  Pj,  GR, 

gain,  receiving  antenna  gain,  transmitter  and  receiver  cable  losses  and  shading 
factor  respectively.  The  shading  factor,  SF,  reflects  the  attenuation  caused 
by  propagation  around  the  fuselage  and  also  includes  a  factor  representing  the 
diffraction  effect  of  a  wing  separating  the  two  antennas.  The  results  are 
based  on  the  work  in  [31]  with  similar  results  appearing  in  [291.  The  last 
term  iti  (4*2)  is  the  Friis  transmission  equation  where  is  the  wavelength  of 
the  chosen  worst  case  frequency  and  D  is  the  shortest  distance  around  the  fuse¬ 
lage  between  the  two  antennas. 

The  selection  of  the  worst  case  frequency  is  performed  in  the  following 
manner.  One  may  expand  (4-1)  and  (4-2)  into 

EH  -  PT  -  20  log1()  f  -  SR  +  K  (4-3) 

where  K  is  a  constant  for  fixed  GL,  G  ,  CL  ,  CL  and  D  and  the  frequency  depend- 
ence  of  SF  is  neglected.  The  transmitter  power,  PT,  is  specified  in  Figure  36a, 
the  selectivity  of  the  receiver  is  specified  in  Figure  36b,  and  the  term  20  log^f 
represents  a  20  dB  per  decade  slope.  The  term  K  in  (4-3)  is  set  to  zero  and  the 
remaining  three  graphs  are  combined  as  in  (4-3)  to  determine  the  worst  case  coup¬ 
ling  frequency. 

The  gain  functions  are  approximated  by  a  two-level  model;  main  beam  gain 
and  mean  side  lobe  gain.  Formulas  are  included  for  computing  these  gains  in 
the  absence  of  known  data. 

The  program  output  consists  of  the  maximum  and  minimum  coincident  fre¬ 
quencies  with  the  selected  worst  case  coincident  frequency.  The  EMI  margin  for 
each  transmitter-receiver  pair  is  printed  out. 

The  field-to-wire  model  is  based  on  the  results  of  [24]  and  translates  an 
incident  electromagnetic  field  into  currents  at  either  end  of  the  line  using 
the  known  line  terminations  and  the  distributed  model  for  the  line  as  described 


GR,  CI^,  CLr,  SF  are  the  transmitter  power,  transmitting  antenna 
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PT  - 
(dBm) 


SR  ^ 
(dBm) 


20log,0f 


in  II.  This  model  is  primarily  intended  for  use  in  predicting  antenna-to-wire 
coupling  through  fuselage  apertures  such  as  gear  doors,  canopies,  etc.  In  this 
useage,  the  antenna-to-antenna  analysis  program  is  used  to  determine  the  power 
incident  upon  the  wire  pair  due  to  an  on  board  transmitter.  This  incident 
power  is  converted  to  electric  field  strength  assuming  plane  wave,  free  space 
propagation  as  described  in  11  and  a  worst  case  field  orientation  is  used  for 
the  f ield-to-wire  coupling  analysis.  Allowances  are  made  for  shielding  of  the 
wires  based  on  empirical  data.  The  data  output  for  this  portion  of  the  program 
consists  of  the  EMI  margin  for  the  receptor  wire  and  each  transmitter.  The 
total  effect  of  all  transmitters  on  the  receptor  wire  is  not  obtained. 

This  program  also  has  a  routine  for  determining  box-to-box  coupling.  This 
routine  only  considers  power  frequency  magnetic  box-to-box  coupling  at  400  H2 
as  well  as  the  third  harmonic,  1200  Hz.  The  transformer  is  considered  as  an 
N  turn  loop  with  current  I.  The  expression  for  magnetic  flux  produced  by  this 
loop  is  easily  obtained  from  elementary  field  theory  (see  [6l).  The  receptor 
is  considered  to  be  a  power  transformer  feeding  the  input  to  a  high  gain  audio 
amplifier.  Consequently,  the  model  of  the  receptor  consists  of  an  M  turn,  open- 
circuited  coil  and  the  induced  voltage  in  this  coil  is  compared  to  the  suscepti¬ 
bility  threshold  of  the  amplifier  in  determining  the  EMI  margin. 

In  this  program,  wire-to-wire,  antenna-to-wire,  antenna-to-antenna  and 
box-to-box  coupling  are  considered  directly.  However,  one  may  also  model  box- 
to-antenna,  box-to-wire,  antenna -to-box  with  a  rewriting  of  the  program.  How¬ 
ever,  there  exists  no  model  capability  to  determine  wire-to-antenna  and  wire- 
to-box  coupling. 

In  using  this  program,  one  must  be  careful  to  determine  the  total,  or  at 
leant  the  "significant,"  energy  incident  upon  a  receptor  from  all  sources.  For 
example,  the  voltage  appearing  across  the  input  terminals  of  a  two  terminal 
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receptor  must  be  the  total  voltage  due  to  coupling  from  all  wires,  cases  and 
antennas  (and  ir,.tersystem  sources  if  significant)  if  one  wishes  to  relate  the 
results  to  compatibility.  If  the  coupling  medium  is  assumed  linear,  then  the 
program  can  still  be  executed  properly  in  modular  form  to  determine  the  energy 
incident  on  a  receptor  due  to  each  source  or  groups  of  sources  and  these  can 
then  be  added  off-line  to  determine  the  total  energy  incident  upon  the  receptor. 
Then  the  results  can  be  related  to  compatibility. 

As  was  discussed  in  III,  the  independent  use  of  the  wire-to-wire  program 
implicitly  assumes  that  the  wire  generators  do  not  interact  with  box  and  antenna 
generators  and  wire  receptors  are  assumed  to  not  interact  with  box  and  antenna 
receptors.  Also  inherent  in  each  of  these  modular  subprograms  (with  the  excep¬ 
tion  of  the  wire-to-wire  routine)  are  the  assumptions  that  the  generators  and 
receptors  other  than  the  pair  being  considered  do  not  interact  with  aach  other 
or  the  pair  under  consideration. 
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V.  THE  TRW  PROGRAM 


We  will  discuss  the  philosophy  of  the  program  first  and  then  consider  the 
models  used.  This  program  represents  a  significantly  new  philosophy  for  the 
EMC  community.  As  was  discussed  in  1,  the  previous  method  of  attack  for  con¬ 
sidering  intrasystem  EMC  consisted  of  applying  somewhat  rigid  and  rather  gen¬ 
eral  specifications  to  the  individual  subsystems  (e.g.,  power  and  signal  lines, 
antennas,  etc.)  comprising  the  total  system.  It  vac-  pointed  out  that  these 
specifications  were  inadequate  in  many  cases.  The  coupling  between  these  sub¬ 
systems  has  been  considered  generally  only  in  an  indirect  way  by  adjusting  the 
specification  limits  through  waivers.  This  adjustment  is  in  most  cases  based 
on  engineering  judgement.  The  difficulties  in  this  regard  are  that  incompati¬ 
bilities  may  still  surface  during  system  test  as  a  result  of  the  human  incapa¬ 
city  to  predict  precisely  the  highly  complex  coupling  between  the  many  sub¬ 
systems.  Also,  the  MIL-STD  specifications  do  not  easily  admit 
trade-offs  in  the  sense  of  required  suppression,  which  would  result  in  mini¬ 
mum  system  weight  and  complexity. 

The  TRW  program  develops  subsystem  specifications  which  are  tailored  to 
the  individual  systems.  Obviously,  each  system  will  have  different  degrees  of 
coupling  between  subsystems  due  to  the  different  physical  configurations  and 
types  of  coupling  paths  between  subsystems.  Then  if  one  determines  analytically 
through  a  mathematical  system  model,  which  predicts  this  coupling,  the  signal 
levels  coupled  between  subsystems,  one  should  be  able  to  arrive  at  specifica¬ 
tions  on  the  subsystems  unique  to  each  system  which  will  achieve  system  com¬ 
patibility  with  a  minimum  amount  of  unneeded  suppression  and  restriction.  The 
TRW  approach  is  directed  along  these  lines. 
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Of  course,  one  may  point  out  that  there  are  many  systems  whose  configura¬ 
tion  and  composition  are  not  fixed.  Also,  many  subsystems  ("black  boxes")  are 
procured  for  general  use  in  different  types  of  systems.  At  first  glance,  this 
may  seem  to  limit  the  program's  applicability  to  fixed  systems  and  subsystems 
designed  for  specific  systems.  Reference  [2*]  contains  valid  arguments  showing 
that  this  need  not  be  necessarily  true.  The  first  argument  against  tailored 
specifications  is  that  particular  types  of  subsystems  are  procurer*  for  use  in 
many  different  types  of  systems  and  these  subsystems  adhere  to  a  fixed  set  of 
specifications  instead  of  a  set  for  each  application.  Tailoring  specifications 
can  still  be  used  here  since  each  particular  system  could  be  modeled  and  the 
composite  minimum  of  the  subsystem  specifications  generated  for  each  system  could 
be  applied  in  the  procurement  of  each  class  of  subsystem  which  would  guarantee 
their  compatibility  with  these  varied  types  of  systems.  On  the  other  hand,  if 
one  does  not  wish  to  model  each  system  and  determine  a  set  of  subsystem  speci¬ 
fications,  then  the  actual  subsystem  test  data  could  be  utilized  to  determine 
if.  these  subsystems  would  be  able  to  integrate  compatibly  with  each  system 
considered. 

The  second  argument  against  tailored  subsystem  specifications  which  are 
unique  to  each  system  is  that  in  many  cases  the  system  environment  is  not  fixed 
or  deacribable,  e.g.,  an  overflying  ai'  raft.  In  this  case,  it  is  impossible 
to  determine  if  any  subsystem  limits  would  insure  compatibility  and  the  best 
"insurance"  that  one  could  obtain  would  be  to  suppress  all  undesired  subsystem 
emissions  and  to  reduce  subsystem  susceptibilities  to  the  state'ut-the-art. 

This  would  be  the  best  we  could  do;  however,  this  is  clearly  not  cost  effective. 
In  this  case,  an  approximation  of  the  environment  would  be  better  than  none  at 
all. 
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Let  us  now  discuss  the  TRW  approach.  A  communication-electronic  system 
consists  of  two  types  of  sources  of  electromagnetic  energy,  desired  or  func¬ 
tional  and  undesired  or  nonfunctional.  An  example  of  a  functional  source  is 
the  fundamental  frequency  of  a  transmitter  whereas  the  harmonics  of  this  funda- 
mental  frequency  are  examples  of  nonfunctional  sources.  If  we  assume  the  coup¬ 
ling  paths  between  sources  and  receptors  are  linear  and  time-invariant,  then  a 
frequency  domain  description  of  these  paths,  their  inputs  and  outputs  is  equi¬ 
valent  to  a  time  domain  description  through  the  Fourier  Transform.  For  illus¬ 
tration,  let  us  suppose  that  the  system  consists  of  p  functional  generators,  q 
functional  receptcrs  and  one  nonfunctional  generator  and  receptor  each,  repre¬ 
sented  by  G^, . . *  »Gp»  Rj . Rq»  «nd  G^,  R^  respectively  as  illustrated  in  Fig¬ 

ure  37.  Let  us  make  the  assumption  that  the  receptors  do  not  "load"  the  gen¬ 
erators.  Also  assume  that  the  generators  do  not  interact  and  the  receptors  do 
not  interact.  With  this  assumption,  we  may  represent  the  coupling  paths  as 
frequency  domain  transfer  functions,  H^Oe).  For  example,  H2N(u>)  represents 
the  transfer  function  between  and  R^. 

The  outputs  of  all  functional  generators  are  frequency  spectrums  repre¬ 
senting  the  Fourier  '’’ransform  of  the  time  domain  outputs  of  these  generators. 
The  phase  information  is  dropped  and  only  the  magnitudes  of  these  Fourier 
Transforms  are  used  yielding  a  worst  case  analysis.  The  susceptibility  of  the 
fur  '-.ional  receptors  is  represented  quantitatively  by  a  single  number  AS^. 

The  susceptibility  index  AS^  is  determined  for  a  particular  receptor  by  assum¬ 
ing  that  the  receptor  will  respond  to  the  sum  of  the  incident  frequency  com¬ 
ponents,  i.e., 


AS  -  f  2  P  (f)  t(f)  dt 

j  fi  ] 


(5-1) 


.10 


where  f  is  frequency,  f^  -  10Hz,  f2  *  10  Hz,  P^(f)  Is  the  "bandwidth"  of  the 
receptor  port  and  f(f)  is  the  total  spectrum  incident  at  the  receptor  port. 
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In  practice,  AS  ^  can  be  more  easily  determined  from  design  considerations  since 
the  ports  for  functional  receptors  have  well  known  properties,  e.g.,  the  input 
co  a  receiver,  as  opcosed  to  nonfunctional  receptor  ports. 

The  program  first  d<  , ermines  functional  compatioility  between  all  functional 
generators  and  each  functional  receptor  through  a  surplus  margin  As'  from 


-  ASj  -  Jfi  pj(f)  {  £  Cl<f)  } 


(5-2) 


If  AS^  <  0  then  the  j-tli  functional  receptor  will  be  rendered  susceptible  due 
to  functional  sources.  Then  design  modifications  must  be  instituted  since  this 
functional  receptor  will  certainly  be  interfered  with  when  all  sources  (function¬ 
al  as  well  as  nonfunctional)  are  in  operation. 

If  all  functional  receptors  are  compatible  with  all  functional  generators 
then  we  will  have  q  surplus  margins  AS^  for  j-l,...,q  which  are  all  greater 
than  or  equal  to  zero.  The  interference  from  the  nonfunctional  source  can 
be  allotted  to  these  surplus  margins  in  the  following  manner.  Compute  for  j*l 


ASX  «  AS'  -  j*  P^f)  HNl(f)  GN<f)  df 
fl 


(5-3) 


If  AS^  is  regative  then  reduce  the  spectrum  of  GN(f}  (which  is  chosen  on  a 
worst  case,  arbitrary  basis)  so  tht»t  (5-3)  will  now  be  nonnegative.  Obviously 

there  are  many  such  reductions  and  spectrums  of  G  (f)  which  will  yield  a  non- 

"  (1) 
negative  AS^  .  Call  thi:  spectrum  G^  (f). 


as2  -  as2  - 


2  -  P2(f)  HN2<f)  GN(f)  df 

1 


(5-4) 


If  AS 2  is  negative  then  reduce  the  spectrum  of  GN ( f )  so  that  AS2  is  nonnegative. 

(2) 

Call  this  resulting  spectrum  G^  (f).  Continue  for  each  j«*l,...,q  to  yield 

spectrums  (f ) , . . .  ,G.^  (f) .  For  each  frequency  (from  a  set  of  selected  fre- 

quency  points)  between  f.  and  f  ,  determine  the  minimum  of  G^^  (f) , . . .  ,G^  (f) , 

1  z  N  N 
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Construct  the  spectrum  GN  (f)  from  these  points.  It  then  is  clear  that  each 

NT 

functional  receptor  will  be  compatible  with  all  of  the  functional  generators 
and  the  single  nonfunctional  generator. 

Now  determine 

V«  -  {  V£>  Vf>} +  <W£)  V£)  <5'5> 

To  measure  compatibility  for  this  system  then,  we  must  tert  the  nonfunctional 
receptor  for  susceptibility  by  applying  R^(f)  to  this  receptor  in  a  laboratory 
situation. 

For  practical  implementation,  the  transfer  functions  H^.(f)  and  H^(f)  for 
1=1,..., p  and  j=l,...,q  must  be  determined.  To  determine  these,  a  "worst  case" 
nonfunctional  generator  and  receptor  is  assumed.  The  physical  location 
of  this  generatoi  and  receptor  is  determined  on  a  worst  case  basis  for  deter¬ 
mining  the  above  transfer  functions. 

The  difficult  off-line  portions  of  the  analysis  concern  the  determination 
of  a  relevant  and  meaningful  susceptibility  index  AS^  and  also  correct  and 
meaningful  worst  case  generator  and  receptor  transfer  functions  above.  These 
two  areas  are  very  crucial  to  the  generation  of  a  meaningful  set  of  specifica¬ 
tions  (tailored  to  the  specific  system). 

The  question  at  this  point  is  how  the  above  analysis  relates  to  developing 
a  set  of  specifications  to  be  applied  to  all  nonfunctional  generators  -he 
same  way  MIL-STD-461  is  used  to  insure  system  compatibility.  Keep  in  mind  that 
in  the  actual  system  test  and  deployment,  each  receptor  (functional  as  well  as 
nonfunctional)  will  he  subject  to  the  affects  of  all  generators  (functional  as 
well  as  nonfunctional).  Here  we  have  developed  a  specification  which  would 
work  tor  a  system  with  only  one  nonfunctional  generator.  Obviously  this  is  a 
highly  unrealistic  situation.  However,  if  we  divide  the  spectrum  G^  (f)  by 


the  number  of  nonfunctional  generators  and  apply  the  resulting  spectrum  to  each 
nonfunctional  generator',  then  we  have  insured  compatibility  if  all  nonfunctional 
generators  are  concentrated  at  the  location  of  G„  and  are  the  same  type  as  G„. 

If  they  are  not,  then  we  still  will  roughly  add  the  effects  of  all  nonfunctional 
generators  or  each  of  the  functional  receptors.  This  is  of  course  an  appr  .i- 
mation  since  the  transfer  functions  between  each  of  these  nonfunctional  gen¬ 
erators  and  each  functional  receptor  will  not  be  equal  to  However, 

we  arc  getting  closer  to  uniquely  tailored  specifications.  The  TRW  program  in 
fact  uses  this1  method  to  develop  the  individual  specifications  for  nonfunctional 


generators. 

The  individual  subsystem  susceptibility  specifications  are  obtained  in  a 
similar  manner  by  dividing  R^f),  given  by  (5  *5‘« ,  by  a  factor  representing  the 
total  number  of  nonfunctional  receptors.  The  resulting  spectrum  is  then  applied 
in  the  laboratory  to  these  nonfunctional  receptors  to  determine  if  compliance 
to  the  specifications  will  relate  to  compatibility. 

It  should  be  pointed  out  that  the  program  is  quite  flexible  since  as  the 
system  design  progresses,  test  ata  on  the  nonfunctional  generators  and  recep¬ 
tors  will  become  available  and  can  be  brought  in  as  data  input  with  the  result 
that  the  program  can  be  used  solely  as  an  interference  anlysis  program  by 
classing  all  generators  and  receptors  (functional  and  nonfunctional)  as 
"functional"  ones.  Also,  once  the  system  has  been  modeled,  then  a  rapid  assess¬ 
ment  of  the  impact  of  granting  waivers  as  well  as  future  system  modifications 
can  be  assessed  with  a  rapid  turn-around  time  since  the  major  portion  of  the 
off-line  work  has  been  done,  i.e.,  the  modeling  of  the  system. 

We  should  point  out  that  one  of  the  difficult  problems  in  implementing 


this  program  will  he 


the  determination  of  receptor  interference  margins. 
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and  input  banJwidths ,  P,v«).  Note  that  we  do  not  determine  these  parameters 
for  nonfunctional  receptors  in  the  early  design  stages.  These  parameters  are 
much  easier  to  determine  for  functional  receptors  than  for  nonfunctional  ones 
in  the  early  design  stages  since  the  functional  receptors  have  a  great  deal 
more  information  about  their  input  ports  and  susceptibilities  than  nonfunctional 
ones  (e.g.,  the  input  to  a  receiver  as  opposed  to  the  power  line  input  for  a 
device) . 

Coupling  Models: 

All  coupling  determinations  are  performed  on  a  generator-receptor  pair 
basis.  For  example,  the  coupling  between  a  wire-type  generator  and  a  wire- 
type  receptor  is  performed  by  solving  only  the  -ircuit  consisting  of  this  pair. 

The  wire-to-wire  coupling  analysis  for  a  generator-receptor  pair  is  deter¬ 
mined  in  two  parts.  The  mutual  capacitive  coupling  between  the  wires  is  deter¬ 
mined  by  neglecting  mutual  inductance.  Next,  the  mutual  inductance  coupling 
is  determined  by  neglecting  the  mutual  capacitance  coirpling.  The  total  coupling 
is  then  considered  to  be  the  sum  of  these  coupling  effects.  Note  that  this  is 
an  approximation  to  the  exact  analysis  of  determining  the  total  coupling  by 
simultaneously  solving  for  inductive  and  capacitive  coupling.  Also,  as  dis¬ 
cussed  above,  only  generator-receptor  pair  calculations  are  made  yielding  first 
order  effects  (i.e.,  the  parasitic  nature  of  other  wires  is  neglected).  In 
this  computation  the  source  data  is  a  magnitude  vs.  frequency  voltage  or  cur¬ 
rent  source  spectrum.  Also,  only  the  magnitude  of  the  transfer  function  is 
used  in  the  coupling  calculation  yielding  a  worst  case  result  for  this  gen¬ 
erator-receptor  pair.  A  significant  modeling  diiference  exists  between  this 
wi.re-to-wire  coupling  computation  and  those  of  III  and  IV  in  that  the  distri¬ 
butive  nature  of  the  generator  and  receptor  wires  are  not  considered  at  all, 
i.e.,  neither  the  distributed  model  nor  the  lumped  h  model  is  used  for  the 
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wires  themselves.  Essentially  this  assumes  that  the  voltages  of  the  load  and 
source  output  of  the  generator  line  are  equal  and  the  currents  of  the  load  and 
source  output  of  the  generator  line  are  equal.  Similar  remarks  apply  to  the 
receptor  line.  At  present  the  load  and  source  impedances  are  considered  to  be 
real  valued  resistances  with  a  capability  to  include  shunt  capacitance.  The 
program  is  not  presently  capable  of  accepting  general  complex  impedances  for 
source  and  load.  Various  shielding,  grounding,  ana  twisting  configurations  are 
included  in  the  capability. 

The  antenna  models  are  simply  the  specification  of  the  E  and  H  fields  at 
the  antenna  site  (or  1  meter  a way)  and  the  replacement  of  the  antenna  with  this 
plane  wave.  The  input  data  needed  is  shown  in  Figure  38.  So  essentially  the 
antenna  is  considered  to  be  an  it  owropic  point  source  radiator  whose  fields  propa¬ 
gate  outward  uniformly  in  all  directions  and  whose  frequency  dependence  is  as  it 

Figure  38.  One  would  have  to  specify  f , ,  f„,  and  E  .  and  E  off-line.  The 
6  r  1  1  2  min  max 

antenna  H  field  is  specified  in  the  same  way  and  both  are  specified  one  meter 
from  the  antenna. 

For  wire  radiation  (for  example,  wire-to-antenna  coupling  and  wire-to-wire 
coupling  at  large  distances)  the  H  fielJ  at  one  merer  from  a  wire  parallel  to 
a  ground  plane  is  determined  in  terms  of  the  current  on  the  wire  using  elemen¬ 
tary  field  theory.  Similarly,  the  E  field  at  one  meter  from  the  wire  is  deter¬ 
mined  in  terms  of  the  wire  voltage. 

The  field  transfer  models  simply  relate  the  above  fields  at  one  meter  from 
the  generator  (antenna  or  wire)  to  the  field  at  the  receptor  (antenna  or  wire) 
as  in  Figure  39.  The  field  transfei  represents  a  decay  factor  determined  off 
line  to  relate  the  field  F^  at^to  the  field  at(l);  for  example, 
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The  field  reception  models  relate  the  fields  at  the  receptor  to  the 
voltage  across  the  load  for  wire  receptors  and  at  the  base  of  the  antenna  for 
antenna  type  receptors.  The  sntenna  reception  essentially  uses  a  frequency 
selection  of  the  form  in  Figure  38. 

In  this  program,  the  wire-to-wire  coupling  is  calculated  on  a  generator- 
receptor  pair  basis  by  adding  capacitive  and  inductive  transfer  for  each  gen¬ 
erator-receptor  pair.  For  wire  pairs  separated  over  moderate  to  large  dis¬ 
tances,  the  E  and  H  field  coupling  is  computed  for  each  generator-receptor  pair. 
Also,  the  antenna-to-antenna ,  antenna-to-wire,  antenna-to-box,  box-to-wire, 
box-to-box,  box-to-antenna,  wire-to-box,  and  wire-to-antenna  coupling  is  cal¬ 
culated  on  a  generator-receptor  pair  basis. 
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VI.  THE  SACHS -FREEMAN  PROGRAM1 

This  work  was  procurred  by  the  Electronic  Sy^ems  .'ivision  (ESD)  of  the 
U.  S.  Air  Force.  We  will  discuss  the  portion  of  the  program  dealing  with  the 
intrasystem  compatibility  investigation  via  compliance  to  MIL-STD-461. 

This  portion  of  the  program  is  focused  on  determining  whether  subsystem 
compliance  to  MIL-STD-461  will  insure  overall  system  compatibility.  The  inves¬ 
tigation  consists  of  three  parts. 

The  Class  I  investigation  determines  box-to-wire  and  wire-to-wire  coupling 
and  uses  the  wire-to-wire  and  field-to-wire  coupling  models  developed  by  TRW 
and  described  in  V.  This  analysis  considers  wire-to-wire  and  box-to-wire  gen¬ 
erator-receptor  pairs. 

For  wire-to-wire  coupling,  the  routine  selects  the  applicable  conducted 
narrowband  limit  (from  CE01,  CE02,  CEO?  or  CE04)  at  the  frequency  of  the  recep¬ 
tor.  This  conducted  limit  represents  a  maximum  allowed  value  of  current  on  the 
emitter  line  and  is  coupled  via  inductive  transfer  to  the  receptor  line  result¬ 
ing  in  a  voltage  across  the  input  to  the  receptor,  vNCj«  This  voltage  is  a 
single  numerical  value  and  not  a  frequency  dependent  spectrum  since  only  the 
receptor  frequency  was  considered.  Next,  the  routine  divides  the  value  of  the 
above  narrow  band  conducted  limit  (obtained  by  evaluating  the  appropriate  CE 
limit  at  the  receptor  frequency)  by  the  emitter  line  impedance  resulting  in  a 
vcltage  on  the  emitter  line.  It  is  not  clear  which  emitter  line  impedance  is 
considered;  for  example,  the  emitter  line  load  impedance  or  the  impedance  seen 
by  the  source.  This  emitter  line  voltage  is  coupled  to  the  receptor  line  through 
capacitive  transfer  resulting  in  a  voltage  across  the  input  to  the  receptor, 

V  .  Again,  V  is  a  single  number  and  not  a  frequency  dependent  spectrum 
NVjVj  NCw 

since  only  the  receptor  frequency  was  considered. 
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The  above  calculations  are  repeated  using  the  applicable  broadband  liaiits 

in  MIL-STD-461  resulting  in  V  and  V  ,  the  inductive  and  capacitive  coupled 

BCI  BCC 

broadband  limits  as  coupled  to  the  receptor  line  respectively.  The  broadband 
limits  are  tnul:iplied  by  the  bandwidth  of  the  receptor  circuit  in  making  these 
calculations. 

The  values  of  the  radiated  emission  limits  (RE01,  RE02)  are  determined 
arparently  at  the  frequency  of  the  receptor.  Each  wire  coupling  pair  separated 
at  moderate  distances  (those  not  considered  above)  and  each  cas_— to-wire  coupling 
pair  are  now  considered  by  transferring  the  values  of  RE01  or  RE02  selected  at 

I 

the  receptor  frequency  to  the  receptor  wire  through  the  TRW  field-to-wire  coup¬ 
ling  models  resulting  in  voltages  at  the  receptor  input  port  V  ,  V  due  to  the 

narrowband  limits  and  V  due  to  the  narrowband  H  field  emission  limit,  RE01. 

BH 

The  above  signal  levels  for  narrowband  and  broadband  interference  are  added 
together  to  produce 


v  ■  V  +  V  +  V 
B  BCI  BCC  BH 

and 


(6-1) 


V-V  +  V  +  v  +v 
N  NCI  NCC  NE  NH 


(6-2) 


where  is  the  broadband  coupled  interference  and  VN  is  the  narrowband  coupled 
interference  each  of  which  is  due  to  wire  and  box  radiators  only  and  are  given 
in  volts.  Each  value  is  a  single  real  number. 

At  this  point,  one  must  specify  a  threshold  voltage  AS  at  which  the  recep¬ 
tor  malfunctions.  This  quantity  is  a  single  number  given  in  volts.  The  sus¬ 
ceptibility  of  each  wire  receptor  to  all  wire  and  box  emitters  only  is  determined 
by  summing  the  total  narrowband  received  voltage  at  the  receptor  input,  V  ,  and 

the  total  broadband  received  voltage  at  the  receptor  input,  VR  .  The  suscepti- 

°T 


bility  of  this  wire  receptor  is  determined  from 
ACTN  »  AS (CM)  -  Vn 


(6-3) 
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and 


ACTB  *  AS (CM)  -  V_ 
®T 


(6-4) 


where  CM  is  a  safety  margin.  If  ACTN  or  ACTB  is  negative,  then  the  wire  recep¬ 
tor  is  said  to  be  susceptible. 

The  second  class  considers  antecna-to-wire  coupling.  The  electric  field 

incident  upon  a  wire  receptor  due  to  a  transmitting  antenna  is  calculated  with 

T39] 

the  EPM-1  path  loss  model  and  is  ’'hen  converted  to  receptor  input  wire 
■oltage  with  the  TRW  field-to-wire  reception  model  as  follows.  The  power  re" 
ceived  by  a  hypothetical  isotropic  antenna  in  free  space  at  the  receptor  wire 
location  is  calculated  through 


P_  -  10  log  P_  +  G_  +  LPROP  +  30  (6-5) 

K  T  T 

where  P  is  the  received  power  of  this  isotropic  antenna,  P  is  the  power  of 
R  T 

the  transmitter  at  its  midband  operating  frequency,  is  the  _,ain  of  the  an¬ 
tenna  in  the  direction  of  the  wire  and  LPROP  is  the  value  of  the  free  space 
path  loss  between  the  antenna  and  the  wire  receptor.  This  received  power  is 
converted  to  electric  field  strength,  E,  using  the  assumption  of  uniform  plane 
wave  propagation  as  described  in  II.  The  TRW  E  field-to-vire  transfer  model 
is  then  used  to  obtain  the  resulting  voltage  across  the  receptor  terminals  due 
to  the  effect  of  this  transmitter.  These  values  are  then  compared  to  the  sus¬ 
ceptibility  threshold  number,  AS,  for  the  wire  receptor.  Apparently  the  effect 
of  each  transmitter  on  the  receptor  is  compared  to  AS  separately  and  if  any  one 
effect  exceeds  AS  then  the  transmitter-receptor  wire  pair  is  said  to  be  incom¬ 
patible  or  at  least  a  potential  problem  area  exists. 

Class  three  represents  wire,  box  and  transmitting  antenna-to-box  coupling. 
For  all  emitter-to-box  combinations,  the  program  calculates  the  E  and  H  field 
signal  levels  at  the  receptor  case  and  tabulates  these.  It  is  not  clear  whether 
this  tabulation  is  the  sum  of  the  effects  due  to  each  emitter  or  whether  the 
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effects  are  tabulated  individually.  Furthermore,  it  is  not  clear  whether  these 
tabulations  are  levels  of  received  signals  versus  frequency  or  at  a  single  fre¬ 
quency;  if  so,  what  frequency  is  used.  No  further  calculations  are  made  and 
apparently  one  would  be  expected  to  perform  laboratory  tests  on  the  box  of 
interest  by  simulating  these  tabulated  field  strength (s). 

Sow  discussion  of  this  program  is  in  order.  In  the  calculation  of  wire- 
to-wire  and  box-to-wire  coupling,  the  value  of  the  applicable  limit  of  MIL-STD-461 
is  selected  at  the  wire  receptor's  frequency.  Zt  is  not  clear  that  all  wire 
type  receptors  are  susceptible  at  only  one  frequency.  For  example,  a  device 
is  not  necessarily  susceptible  to  signals  entering  the  power  line  input  at  only 
the  power  supply  frequency.  To  perform  a  proper  analysis  m  this  regard,  one 
should  determine  the  coupled  signal  levels  at  either  all  frequencies  or  only 
those  for  which  the  device  is  susceptible  at  the  wire  input  port  (if  such  a 
determination  can  be  made). 

Finally,  the  modular  useage  of  this  program  will  point  out  only  certain 
areas  of  susceptibility.  For  example,  if  one  determines  a  wire  type  receptor's 
port  voltage  due  to  all  box  and  all  other  wire  couplings  and  if  this  vah.e  ex¬ 
ceeds  AS  (the  susceptibility  value  of  the  receptor)  then  one  can  surely  state 
that  this  receptr;  will  be  interfered  with.  This  is,  of  course,  completely 
dependent  upon  a  meaningful  susceptibility  value,  AS.  However,  if  this  value 
does  not  exceed  AS,  then  can  one  say  that  this  wire  receptor  will  operate  com¬ 
patibly  with  the  system?  Obviously  this  statement  cannot  be  made  with  confi¬ 
dence.  This  is  due  to  two  considerations.  First  of  all,  this  receptor  may  be 
susceptible  to  signals  other  than  the  single  frequency  chosen  for  the  analysis. 

We  have  not  considered  incident  signals  at  other  frequencies.  The  second  reason 
is  a  result  of  the  fact  that  we  have  not  added  in  the  antenna  and  box  coupled 
effects.  Even  if  there  exist  no  antennas  within  the  system  whose  fundamental 
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frequency  is  the  same  as  the  chosen  wire  receptor  frequency,  incompatibilities 


may  result  since  this  wire  receptor  may  be  susceptible  to  frequencies  other 
than  the  one  designated  as  the  "receptor  frequency"  and  in  particuJ  ,r  one  of 
the  transmitter  fundamental  frequencies.  Also,  one  should  consider  not  only 
the  fundamental  frequency  of  the  transmitter  but  also  the  spurious  responses. 


VII.  CONCUTTON 


In  conclusion,  we  may  review  our  discussions  with  a  brief  comparison  of 
the  philosophies  and  methods  of  attack  in  each  program. 

First  of  all,  we  might  point  out  that  none  of  the  above  programs  consider 
intersystem  effects.  Of  course,  the  modeling  and  inclusion  of  intersystem 
effects  can  be  a  difficult  problem,  yet  they  will  effect  each  system  in  an 
actual  deployment  and  should  be  considered. 

As  for  philosophies,  basically  the  Boeing,  McDonnell-Douglas,  and  Sachs- 
Freeman  programs  consist  of  an  analytical  capability  to  determine  the  electro¬ 
magnetic  energy  coupled  between  subsystems.  This  capability  is  used  to  deter¬ 
mine  whether  or  not  compliance  to  the  subsystem  specifications,  e.g.,  MIL-STD-461, 
will  ensure  system  compatibility.  This  determination  is  crucially  dependent 
upon  having  a  meaningful  quantitative  measure  of  the  susceptibility  of  a  device. 
The  Boeing  and  McDonnell-Douglas  programs  use  a  frequency  dependent  suscepti¬ 
bility  threshold  as  in  the  MIL-STD's  which,  as  was  pointed  out  in  the  intro¬ 
duction,  needs  justification  since  devices  basically  respond  to  total  incident 
energy.  The  use  of  a  broadband  susceptibility  criterion  provides  only  a  rough 
measure  of  this.  The  Sachs-Freeman  program  appears  to  consider  a  device  sus¬ 
ceptible  at  only  one  frequency  which  seems  rather  unrealistic.  The  TRW  program 
on  the  other  hand  considers  a  device  susceptible  to  a  measure  of  total  incident 
energy  since  the  effect  of  an  incident  spectrum  is  summed. 

As  for  models  of  the  coupling  paths  between  subsystems,  all  programs  are 
in  universal  agreement  as  to  the  use  of  a  frequency  domain  description  of  the 
coupling  paths  as  opposed  to  a  time  domain  description.  As  discussed  in  the 
introduction,  this  seems  to  be  a  valid  approach.  The  wire  coupling  models  use 
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lumped  approximations  to  the  exact  distributive  calculation.  Apparently  this 
approximation  is  considered  to  outweigh  the  suspected  increase  in  computation 
time  with  the  more  exact  distributive  model.  These  approximations,  however, 
only  are  valid  for  frequencies  where  the  circuit  dimensions  are  much  less  than 
a  wavelength  and  it  has  not  been  conclusively  demonstrated  that  the  distributive 
calculation  necessitates  the  large  increase  in  computational  i  ’me  over  the 
lumped  approximations  as  seems  to  be  expected. 

All  programs  essentially  perform  the  calculations  of  transferred  energy 
on  an  emitter-rec-tptor  pair  basis  and  assume  that  other  generators  and  receptors 
will  not  affect  the  pair  being  considered  in  the  calculations.  This  obviously 
reduces  the  required  computation  time  substantially  since  the  size  of  the  inves¬ 
tigated  circuits  (for  wire-to-wire  coupling)  generally  requires  the  minimum 
computational  effort.  The  one  exception  to  this  generator-receptor  pair  cal¬ 
culation  is  the  McDonnell-Douglas  wire-to-wire  coupling  program. 

The  antenna  coupling  calculations  also  assume  no  interaction  between 
emitters  and  between  receptors.  For  example,  the  mutual  coupling  between  trans¬ 
mitting  antennas  is  not  considered.  Also  no  computations  of  intermodulation, 
crossmodulation  or  spurious  responses  are  included  for  receivers  attached  to 
antennas.  The  propagation  models  use  the  free  space  path  loss  equations  de¬ 
scribed  in  II  which  rely  heavily  upon  obtaining  accurate  antenna  gain  factors 
(in  the  direction  of  transmission).  Furthermore,  these  gain  factors  are  gen¬ 
erally  assumed  i.  dependent  of  frequency  which  seems  to  be  an  unrealistic 
assumption. 

Box  -adiation,  reception  and  coupling  has  been  minimally  considered  appar¬ 
ently  for  the  reasons  given  in  II. 

the  most  important  point  to  be  brought  to  the  attention  of  users  of  these 
types  ot  piogiams  i.s  the  tol  lowing.  In  the  actual  system  deployment,  each 
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subsystem  will  be  effected  to  some  degree  by  all  sources  (wire,  case,  and 
antenna  coupled  effects  as  well  as  intersyatem  effects).  If  one  determines 
the  effects  on  a  receptor  due  to  wire  type  sources  only  and  obtains  the  result 
that  the  receptor  will  be  interfered  with,  then  it  is  valid  to  conclude  that 
this  receptor  will  be  rendered  susceptible  in  a  system  test.  If,  however,  the 
receptor  is  found  to  be  not  susceptible  to  only  wire  sources,  then  one  cannot 
state  that  the  receptor  will  not  be  interfered  with  during  system  test  si  ice 
the  effects  due  to  othor  sources  (intersystem,  case,  and  antenna)  have  not 
been  added  in.  Considering  partial  effects  of  sources  on  a  receptor  will  not 
necessarily  yield  the  complete  system  compatibility  picture! 
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